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sea  dates  back  to  this  patent  by  C.  Chilowski.  submitted  in  1924.  .Note  the 
sophisticated  4  beam  Janus  design  which  is  still  used  for  first  order  pitch  and  roll 
compensation. 

Fig.  2.2  :  Range  Gated  Shipboard  Doppler  System.  Four  beams  (fore.  aft.  8 

port  arid  starboard)  equally  inclined  from  the  ship’s  vertical  axis  acoustically 

probe  the  ocean.  The  Doppler  shift  in  backscattered  signal  measures  the  relative 
velocity  between  ship  and  ocean.  Range  gating  allows  depth  profiling  of 
currents. 

Fig.  2.3  :  A  Doppler  current  profiler  which  has  been  bottom  mounted  for  11 

vertical  profiling  of  local  currents. 
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Fig.  3.1  :  Depth  range  of  shipboard  measurements  is  shown  by  the  percent  of  17 

pings  for  which  the  profiler  electronics  indicated  valid  returns  at  each  range  bin. 

Two  12  hour  subsets  of  data  are  shown,  one  each  from  1981  and  1982.  Details  of 
such  validity  profiles  vary  with  many  factors  including  location  and  time  of  day. 

Fig.  3.2  :  Variability  in  acoustic  measurements  of  relative  velocity,  shown  as  19 

histograms  of  the  deviations  of  single  ping  measurements  from  100  ping  average. 

Fig.  3.3  :  Autospectra  of  ping-to-ping  relative  velocity  measurements  for  3  21 

subsets  of  data  taken  at  nearly  constant  shipspeed  while  (a)  steaming  into  seas. 
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velocity  F.  dashed  line  is  portward  velocity  P.  Note  peak  at  surface  wave 
encounter  frequency  rising  above  a  white  noise  background. 

Fig.  3.4  :  Coherence  of  F  (solid)  and  P  (dashed)  with  pitch  <t>  and  roll  p  of  the  22 

ship.  Data  in  left  panels  taken  while  on  station  (same  as  middle  panel.  Fig.  3.3); 
data  in  right  panels  taken  while  steaming  into  seas  (  same  as  upper  panel. 

Fig.  3.3) 


Fig.  3.5  Expected  residual  wave  induced  and  white  noise  in  acoustic  23 

measurements  of  relative  velocity  after  block  averaging,  as  a  function  of 
averaging  length. 

Fig.  4.1  :  Sample  time  series  of  measurements.  Top  panel  shows  the  35 

acoustically  determined  relative  velocity  V.  middle  panel  shows  dx0  dt  from 
LORAN-C.  Bottom  panel  shows  the  sum  of  these,  the  fix-to-fix  current 
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Fig.  4.2  :  Spectrum  of  fix-to-fix  current  measurements  for  v  (  solid  line.  36 

alongshore  current.  317"T  )  and  u  (  dashed  line,  cross-shore  current.  47°T  ). 

Fig.  4.3  :  Same  as  Fig.  4.2  except  normalized  by  sin'trfAt.  White  noise  in  36 

LORAN  fixes  yields  a  flat  spectrum  for  fix-to-fix  currents  under  this 
normalizat  ion. 

Fig.  4.4  :  The  mean  square  error  in  filtered  fix-to-fix  current  estimates  43 

(cm*  sec  ")  is  predicted  as  a  function  of  filter  length  (fix  intervals)  for  3  filters 
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in  v  and  divergence  in  u  across  the  shelf. 
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Fig.  6.1  :  Alongshore  component  of  surface  wind  stress  (dynes  cm*)  at  NDBO  61 

46013.  computed  from  hourly  wind  measurements.  The  hourly  stress  values  were 
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50  and  100  fathom  isobaths,  and  are  numbered  consecutively  from  northern 
station.  See  Appendix  for  times  occupied. 
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over  the  observations  of  Fig.  6.3a. 

Fig.  6.3c  :  Standard  deviation  of  fluctuations  in  u(x.z).  v(x.z)  (cm/sec)  for  66 

observations  in  Fig.  6.3a. 

Fig.  6.3d  :  RMS  expected  error  of  the  average  of  u(x.z).  v(x.z)  (cm' sec).  67 

assuming  observations  in  Fig.  6.3a  are  independent. 

Fig.  6.3e  :  Vertical  shear  irt  the  mean  fields  9u/'9z,  9v/9z.  First  differences  of  68 

the  mean  fields  (Fig.  6.3b)  were  smoothed  using  a  (3*  3)  triangular  weight  filter. 

Fig.  6.3f  :  Cross-shore  divergence  Bu/9x  and  shear  9v/9x  in  mean  currents.  69 

calculated  as  in  Fig.  6.3e.  Values  enclosed  in  boxes  are  significantly  different 

from  zero  at  the  95ca  confidence  level,  assuming  errors  in  the  mean  are  given  by 
Fig.  6.3d  and  that  these  errors  are  independent  between  adjacent  bins.  The 
latter  is  almost  certainly  optimistic. 


Fig.  6.4a  :  Average  u(x.z).  v(x.z)  along  each  CTD  line.  Note  the  jet  in  the 
average  alongshore  current  v.  the  core  of  which  intensifies  and  moves  offshore 
from  north  to  south. 


Fig.  6.4b  :  Standard  error  of  the  mean  for  Fig.  6.4a. 


Fig.  6.4c  :  Standard  deviations  of  u.  v  about  the  means  in  Fig.  6.4a. 


Fig.  6.5  :  Alongshore  divergence  ov  dy  between  pairs  of  primary  CTD  lines. 
Values  which  were  non-zero  at  the  95^  confidence  level  as  discussed  in  Fig.  6.3f 
are  boxed. 


Fig.  6.6  :  Principal  axes  of  fluctuations  about  the  mean  currents  on  each  CTD 
station.  Where  fluctuations  are  nearly  equal,  the  orientation  of  the  principal  axes 
is  random. 


Fig.  7.1  :  DAL  currents  at  depth  of  28m.  obtained  during  a  survey  from 
26  April  0045  to  28  April  1400.  1981  (all  times  are  UT).  plotted  over  a  NOAA6 
IR  image  from  27  April  1611.  Tic  marks  every  30  minutes  of  latitude  and 
longitude  (56  and  44  km  respectively)  and  a  scale  vector  for  the  currents  are 
shown.  Pt.  Arena  (38°  57'  N.  123°  44'  W)  serves  as  a  landmark.  Lighter  shades 
correspond  to  colder  water.  For  clarity  of  presentation  all  measurements  within 
a  3  km  radius  have  been  averaged.  Note  the  concentrated  zones  of  flow  in  which 
cold  water  upwelled  at  the  coast  it  carried  offshore.  Such  flow  structures  may 
make  significant  contributions  to  the  property  transport  budgets  in  the  coastal 
zone. 

Fig.  7.2  :  Same  as  Fig.  7.1.  except  for  ship  survey  of  2  May  1800  to  6 
May  0330.  1981  and  NOAA6  image  from  6  May  0345.  Survey  was  made  during 
a  period  of  sustained  strong  equatorward  winds  (Fig.  6.1).  Compare  with  Davis' 
Fig.  5  (1984a). 


Tig.  7.3  :  Ship  survey  from  18  May  1600  to  22  May  0300.  1981.  NOAA6 
image  from  20  May  1550. 

Fig.  7.4  :  Ship  survey  from  23  May  0100  to  25  May  0600.  1981.  NOAA6 
image  from  20  May  1550  (same  as  Fig.  7.3). 


Fig.  7.5  :  Ship  survey  from  20  April  1600  to  22  April  1200.  1982.  NOAA7 
image  from  24  April  2210  (note  delay  from  survey  time).  This  flow  occurred 
during  a  cessation  of  winds  following  a  period  of  strong  southward  wind  forcing, 
and  is  discussed  in  Seel  ion  7.3. 
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Fig.  7.6  :  Ship  survey  from  22  April  1400  to  24  April  1130.  1982.  NOAA7 
image  from  24  April  2210  (same  as  Fig  7.5). 

Fig.  7.7  :  Ship  survey  from  9  July  0300  to  10  July  1730.  1982.  NOAA7  image  89 

from  9  July  2211. 

Fig.  7.8  :  Ship  survey  from  11  July  1100  to  13  July  1400.  1982.  NOAA7  90 

image  from  9  July  2211  (same  as  Fig.  7.7). 

Fig.  7.9  :  Ship  survey  from  14  July  1400  to  16  July  0920.  1982.  NOAA7  91 

image  from  14  July  2254.  Wind  forcing  was  very  strong  and  increasing  during 

the  survey  period.  Note  the  strong  current  magnitudes  and  shears  across  the 
northern  transects. 

Fig.  7.10  :  Ship  survey  from  16  July  1730  to  18  July  0845.  1982.  NOAA7  92 

image  from  15  July  2241.  Survey  was  performed  from  south  to  north:  winds  were 
strong  to  moderate  along  the  Central  and  .North  lines,  but  had  calmed  to  airs 
during  the  surveys  of  the  Arena  and  Irish  Gulch  lines.  Northward  flow  at  the 
inshore  ends  of  the  CTD  lines  is  seen,  although  the  flow  at  the  offshore  end  of 
the  two  northern  lines  remains  strong  and  strongly  sheared.  Note  the  developing 
eddy  along  the  North  line. 

Fig.  7.11  :  Ship  survey  from  19  July  1900  to  22  July  0400.  1982.  NOAA7  93 

image  from  22  July  2259.  Winds  were  calm  during  the  survey  of  southern  3  lines, 
but  were  moderate  to  strong  along  the  northern  3  lines.  The  eddy  along  the 
North  line  has  developed  strongly.  Davis  (1983)  shows  tracks  of  surface  drifters 
which  circulated  in  this  eddy  for  a  week  and  more,  making  several  circuits. 

Figs.  7.12  :  Cross-shelf  sections  of  alongshore  current  with  depth  during  the  94 

wind  relaxation  of  April  1982.  Sections  are  from  the  Five  primary  CTD  lines, 
plus  a  CN  line  midway  between  the  Central  and  North  lines.  Plusses  show 
location  of  measurements.  Contour  interval  is  10  cm  sec.  Equatorward  currents 
are  shaded. 

Fig.  7.13  :  Alongshore  current  v(x)  at  28m  for  the  4  Central  line  occupations  95 

from  20-24  April. 

Fig.  7.14  :  Cross-shelf  sections  of  e'ongshore  currents  with  depth  for  the  four 
surveys  of  the  Central  line  from  20-24  April. 


96 


Fig.  7.15  Time  series  of  moored  measurements  showing  alongshore 
component  of  wind  stress  r  at  NDBO  46013.  and  temperature  difference  and 
dv  dz  between  instruments  at  5  and  10m  depth  at  C3  vs  day  of  1982. 

Fig.  7.16  :  At  least  4  large  tongues  of  cold  water  are  seen  in  this  NOAA6  IR 
image  of  sea  surface  temperature  between  Cape  Mendocino  and  Pt.  Conception 
(image  provided  by  Pierre  Flament). 

Fig.  7.17  :  Currents  at  28m  depth  from  offshore  survey  of  4-10  July  1981. 
overplotted  on  NOAA6  IR  image  of  sea  surface  temperature  from  8  July. 


Fig.  7.18  :  Currents  at  28m  depth  from  offshore  survey  of  22-26  July  1982. 
overplotted  on  NOAA7  IR  image  of  sea  surface  temperature  from  22  July. 

Fig.  7.19  :  Ship  track  during  1982  offshore  survey,  indicating  locations  of 
sections  through  the  cold  water  tongue.  Clouds  obscured  the  area  during  the 
actual  survey  time,  and  the  location  of  the  front  and  detailed  structure  of  the 
temperature  field  were  not  available  from  satellite  IR. 

Fig.  7.20  :  T-S  diagram  from  underway  thermosalinograph  during  jet  survey 
(data  kindly  made  available  by  Pierre  Flament.  Libe  Washburn  and  Larry  Armi). 
Sampling  depth  is  approximately  am. 

Figs  7.21  :  Across-jet  sections  of  along-jet  currents  through  the  cold  tongue 
and  adjacent  waters.  Also  shown  are  the  temperature,  salinity  and  density  (as 
o  f,)  at  5m  depth  from  the  underway  thermosalinograph  (data  kindly  made 
available  by  Pierre  Flament.  Libe  Washburn  and  Larry  Armi). 

Fig.  A.l  :  Central  Line  occupations  vs.  day  of  year. 


Fig.  A. 2  :  Elk  Line  occupations  vs.  day  of  year. 

Fig.  A. 3  :  Irish  Gulch  Line  occupations  vs.  day  of  year. 
Fig.  A. 4  .  Arena  Line  occupations  vs  day  of  year. 


Fig.  A. 5  :  North  Line  occupations  vs.  day  of  year. 


Fig.  A. 6  :  Ross  Lint-  occupations  vs.  day  of  year. 


114 


Fig.  A.T  ;  Fifty  Fathom  Line  occupations  vs.  day  of  year. 


Ill 


Fig.  A.t<  :  Hundred  Fathom  Line  occupations  vs.  day  of  year. 
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[)  \1.  operating  parameter 


1  abie  1  :  Nomina 
transmit  irequer,c\  ;  .'{(id  kHz 


u  a\  ru'ngt  r. 

f. 

l!.f>  cm 

repeal  interval 

(>.63  sec 

pulse  ierigti, 

T 

20  msec 

range  bin 

A' 

10  msec 

beam  angle 

Or, 

30°  from  vertical 

beam  width 

2C  to  3°  from  centerline 

transducer  depth 

Z 

5  m 

Table  3.2  Range  Bin  Geometry 


n  range  bin  number  1.2 32 

e  speed  of  sound 


range  to  head  of  pulse  after  time  t 

c(T-i) 

range  to  tail  of  pulse  after  time  t 

ct 

duration  of  range  bin  ri 

(n-  1) Ar  to  nAr 

For  the  nl  range  bin 

Nominal  value  (m) 

center  range 

(C  2)(-l-(n-i.)Ar) 

3.8  -  7.5n 

center  depth 

(c  -  2)(-21  -  {n - -i-)Ar  )cosoCl-zc 

8.3  —  6.5n 

range  extent  of  returns 

(c  '2j(T*  Ar ) 

22.5 

range  extent  continuous!;,  ensonified 

(<•  2)  Ar 

7.5 

depth  extent  continuously  ensonified 

(c.'2)ArcosQ0 

6.5 

greatest  range  sampled 

(c/2)(T-(n-l)Ar) 

7 . 5  ( n  —  1 ) 

Chapter  3 

DATA  PROCESSING  AND  ERHOK  ANALYSIS 

Making  shipboard  estimates  of  current  profiles  requires  knowing 
(I)  x  It),  the  instrument  position  as  a  function  of  time,  and 
'2l  V(z.t  |  the  relative  velocity 

In  this  section  we  discuss  'me  measurements  made  during  CODE  to  obtain  these  data, 
devoting  special  attention  to  identifying  potential  sources  of  error.  In  the  next  chapter, 
these  data  will  be  used  to  estimate  the  field  of  currents  u ( z .  1 : x , , )  '  u.  Following  that, 
the  currents  inferred  from  the  shipboard  system  will  be  compared  with  moored  current 
meter  measurement"  in  Chapter  5. 

3.1.  Instrumentation 

The  (lata  acquisition  system  for  this  study  was  designed  and  implemented  In 
Lloyd  Kegier.  Tin  acoustic  measurements  were  made  from  the  K  \  Wecoma  using  a 
prototype  of  the  Ametek-St raza  DCP4015.  controlled  from  a  Commodore  I’Ll  21)01 
1 1 or i it  computer.  Operating  characteristics  of  the  acoustic  system  are  given  in  Table  3.1. 
Once  every  0.(>3  set  .  a  300  kll/  acoustic  pulse  of  20  msec  duration  was  transmitted.  The 
Doppler  shift  Hi  tlit'  acoustic  returns  was  measured  by  the  Ametek  electronics  for  each 
beam  in  32  range  bins,  each  of  duration  As  10  msec.  Table  3.2  give-  tin  nominal 
correspondence  between  range  bin  and  physical  depths.  The  measurement'-  in  each 
range  bin  also  contained  a  "quality"  flag  to  indicate  whether  tin-  frequency  detection  loop 
was  locked.  The  ship's  heading  (from  the  ship’s  Sperry  \1k.37  master  gyrocompass), 
ocean  temperature  at  the  rarisducer.  and  time  (from  tin  internal  PET  clock)  were 
recorded  for  each  pulse.  For  every  100  pulse"  (approximately  once  per  minute),  a 
LORAN-(  fix  was  taken  from  a  North-lar  3000  receiver.  In  addition,  during  CODE  1 
(I9M|.  gyroscopic  measurements  of  pit  r  r  i  and  ml!  wen  sampled  1  (>  times  per  ping,  with 
10  msec  between  samples.  \  1 1  data  wen  recorded  without  further  processing  onto  9 
track  magnetic  tape 

3.2.  LORA  N-C 

Measurement"  of  the  "hip"  latitude  A  and  longitude  L  from  a  Northstar  5000 
I.OKAVC  receiver  were  recorded  throughout  the  experiment  at  intervals  of  approxi¬ 
mately  70  seconds.  FOR  \N-('  F  t  lee-  l,<>ng  Range  Aid  to  Navigation  network  main¬ 
tained  by  t  lie  1  ,  S  C  oast  Cuard.  It  employs  low  frequency  (100  kHz)  radio  signals  to 
provide  a  hyperbolic  navigation  grid.  Position  on  this  grid  is  defined  by  the  differences 
in  reception  time-  of  synchronized  signals  transmitted  from  at  least  three  well  separated 
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Pettigrew  and  Irish  (19831  have  made  current  measurements  with  a  four  beam 
acoustic-  Doppler  instrument  similar  in  most  respects  to  the  one  described  above,  the 
major  difference  being  that  their  system  is  designed  to  make  measurements  from  a  sta- 
tionar>  platform.  The  transducer  is  attached  in  an  upward-looking  configuration  to  a 
frame-  which  rests  on  the  sea  floor  (see  Fig.  2.3).  In  Fig.  2.4.  hour!)  averages  of  current 
measurements  from  their  Doppler  system  art-  plotted  along  with  simultaneous  hourly 
average  measurements  from  a  string  of  moored  current  meters  located  approximately 
3(H)m  away  The  striking  agreement  between  the  time  series  indicates  that  Doppler  back- 
scatter  provides  reliable-  measurements  of  ocean  flow  velocities  relative  to  the  transducer 
for  flow  speeds  up  t o  eirder  30  cm  sec. 

Mounting  the  instrument  on  a  movable  platform  such  as  a  ship  introduces  a 
number  of  complications.  In  particular,  to  obtain  ocean  currents  the-  platform  rotation 
and  translation  must  be  removed  from  the  measurements.  This  requires  introducing 
imperfect  measurements  of  the  platform  motion  which  add  error  to  the  final  result. 
Another  important  difference*  is  that  the-  acoustic  environment  at  the-  ship  s  hull  is  much 
metre  variable  than  that  found  at  the  sea  fleior.  Furthermore,  since  the-  ship  cart  move  at 
speeds  which  are  one*  or  two  eirders  of  magnitude  larger  than  the  currents,  signal-to-noise 
limitations  can  be-  severe. 

This  discussiem  has  glossed  over  a  number  of  important  details.  For  example. 
Le)  2.1  is  only  an  approximation  to  the  Doppler  shift:  refractive  effects  from  sound-speed 
gradients  and  current  shears  have  been  neglected  completely,  as  have  second  order  and 
higher  terms  in  V  c.  Numerical  solutions  which  include  these  effects  have  been 
explored  on  a  case  by  case  basis  in  work  with  Regier.  In  general,  neglect  of  these  factors 
does  not  cause  significant  errors  in  the  determination  of  V  when  a  four-beam  Janus  sy  s¬ 
tem  is  used.  The  same  cannot  be  said,  however,  for  other  beam  geometries  such  as  the 
three-beam  Janus  system.  Moreover,  the  effects  of  finite  beam  width,  side  lobe  levels 
and  multiple  scattering  have  been  ignored.  Rather  than  exhaustively  catalogue  the 
potential  errors,  we  shall  turn  to  examination  of  the  measurements  themselves. 
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i"j  -  siiio.-jc  coso,.z  rp  smo(,y  coscipZ  ( 2 .  u ) 

r\  -  sinac.x  -coso0z*  r4  sino0y  coso(,z ' 

The  spatial  dependence  of  the  relative  velocity  V(r)  in  Eq  2.2  comes  entirely  from  shears 
in  the  current  field,  as  mat  be  seen  from  Eq  1.2.  Since  the  currents  are  expected  to  be 
vertical!}  sheared,  we  expand  V(r)  about  its  value  direct!}  below  the  ship.  With 

r  \  x  -  \  y  -  z  £ 

V(r)  V(z')  ■  (r  z  £')'V  V(z')  -  ...  (2.6) 

the  Doppler  shift  (2.2)  at  time  r  can  be  written  for  each  beam  as 


(sino0x  '  -  cosof.z  ' ) 
(sino0y '  -cosa0z  ) 


Af.(-) 

2fp 

— -  v(z'; 

dV(z') 

|  -  z  tartar - ■ — *- 

c 

dx  ' 

Afo(r) 

2fo  v, 

k  dV(z') 

-  Z  tailftr, - 

c 

dy  ' 

Af 3(r) 

2fn 

— -  V(z  ' 

,  dV(z') 

)  -  z  tanop. - - 

c 

dx  ' 

Af4(r ) 

2ft, 

— -  V (z  ' 

.  dV(z') 

z  tano0 — ~ 

c 

dy 

[-stna0x  -coscipZ 


•  ( -  sina0y  '  -  cosa  p£ ' ) 


where  z'(r)  —  (—  -r)cosoc,  is  the  quasi-vertical  coordinate.  If  we  denote  the  com¬ 

ponents  of  V  by  (F.P.H)  as  in  Fig.  2.2.  then 

r,  M  c  Af|(r )- AfjM  ,  dH  , 

2ftl  ismop  dx 

p(/'i  Af<(r)  -°( z' j^)  -•••  <2-8) 

21^1  2siiio  d\ 

r  Af,(r)- Af2(r)- Af,(r)- Af„(r)  .  ,  dF  dP 

H(z  )  — - - - 0(  z  tan  a0(- — ; — -  r)  ) 

2f,  4coso0  dx  dy 


Thus,  by  combining  the  measurements  from  complementary  beams  the  three  components 
of  relative  velocitv  can  be  determined  as  a  function  of  z  '.  the  quasi-vertical  coordinate. 
Turbulence  and  high  frequency  internal  waves  act  as  noise  to  this  determination  by 
adding  short  scale  variability  through  the  terms  on  the  right  hand  side  of  Eq  2.8. 


ACOUSTIC  BEAMS 


Fi  ,22;  Range  Gated  Shipboard  Doppler  System.  Four  lwn»  (fore, ft 

~  TrStUy  between  sh,p 
and  ocean.  Range  gating  allows  depth  profiling  of  currents. 


I 


I 


Then*  is  a  wide  variety  of  acoustic  targets  in  the  ocean  including  zooplankton, 
air  bubbles,  fish  and  the  ocean  bottom  itself.  Hath  such  target  can  contribute  energy  to 
the  backseat  tered  signal  reaching  the  transducer,  at  a  Doppler  shifted  frequency 
corresponding  to  its  own  relative  velocity  V,.  Thus,  in  addition  to  the  uncertainty  in 
frequence  measurement  introduced  by  short  resolution  time  at  each  depth,  the  signal 
heard  at  the  transducer  will  contain  a  spectrum  of  frequencies  due  to  the  range  of  target 
velocities.  If  the  velocitx  of  the  targets  through  the  water  is  random,  then  averaging 
over  many  scattering  events  yields  a  spectrum  which  is  centered  at  a  Doppler  shifted  fre¬ 
quency  corresponding  to  V.  the  relative  velocity  of  the  water. 

Because  the  acoustic  energy  is  transmitted  as  a  short  pulse  of  duration  T.  the 

returns  heard  at  any  given  time  r  after  transmission  must  have  scattered  within  the  par¬ 

cel  of  water  at  ranges  cr  2  to  c(T  -r)  2  along  the  beam.  Returns  heard  at  longer  r  were 
scattered  from  further  down  the  beam.  Thus  by  measuring  the  Doppler  shift  as  a  func¬ 
tion  of  r  it  is  possible  to  profile  the  along-beam  component  of  V  as  a  function  of  along- 
bearrt  distance. 

Af,(r)  —  (V(ri)T*i)  (2.2a) 

c 

r  ^-(  -I  -  r  )r,  (2.2b) 

In  practice,  the  returning  signal  is  range  gated  -  analyzed  in  discrete  blocks  of  time  Ar 
called  range  bins  -  rather  than  continuously  recorded.  This  allows  a  finite  time  to 
resolve  the  Doppler  shift.  Note  that  allowing  longer  Ar  increases  the  accuracy  with 
which  the  spectrum  can  be  resolved,  but  decreases  the  range  resolution. 

The  full  relative  velocity  vector  V  can  be  determined  as  a  function  of  the  verti¬ 
cal  coordinate  by  using  several  beams,  provided  that  V  changes  primarily  along  the  vert¬ 
ical  symmetry  axis  of  the  beams.  This  may  be  seen  by  adopting  a  coordinate  system 
tied  to  the  ship,  with  origin  at  the  acoustic  transducer  and  coordinate  axes  which  point 
along  the  fore,  port  and  heave  directions  x'.y'  and  i  respectively  (see  Fig.  2.2).  In 
these  coordinates  the-  iir  beam  direction  is 

r.  cosfjSineijX  -  siru’jsinojy '  rosojZ '  (2.3) 

where  v,  is  the  azimuthal  angle,  and  o,  is  the  elevation  angle,  of  the  ith  beam.  For  the 
four-beam  Janus  system,  these  angles  are 

o;  o(l  i  1. 2.3.4  (2.4) 

V,  (i  Dy 


so  that 


Chapter  2 


DOPPLER  ACOUSTIC  LOG 

For  t his  study  V(z).  the  relative  velocity  of  the  water  as  a  function  of  depth, 
was  determined  acoustically  using  a  commercially  available  instrument  system  (Rowe 
and  Young.  1979)  which  we  shall  call  the  Doppler  Acoustic  Log.  or  simply  DAL.  The 
remote  measurement  of  velocity  using  a  Doppler  shift  has  a  long  history:  applications  are 
found  in  fields  ranging  from  astrophysics  to  law  enforcement.  Meteorologists  for  some 
time  have  used  the  technique  to  remotely  measure  winds  from  ground  based  radar  sys¬ 
tems  (Lhermitte.  1973).  Since  1972.  R.  Pinkcl  of  S10  has  applied  it  to  make  oceano¬ 
graphic  measurements  from  the  quasi-stationary  research  platform  FLIP.  Most  recently 
Regier  (1982).  Joyce  et  al.  (1982).  and  Joyce  and  Stalcup  (1984)  have  used  the  technique 
to  make  shipboard  current  measurements. 

It  was  realized  as  early  as  1924  that  the  Doppler  shift  could  be  exploited  to 
make  a  "speed  log"  for  ships.  By  measuring  acoustic  backscatter  from  the  ocean  bottom 
or  from  the  water  column  itself  and  determining  the  Doppler  shift  between  transmitted 
and  received  signals,  the  velocity  of  the  ship  relative  to  the  scattering  medium  can  be 
inferred  in  a  manner  to  be  outlined  below.  The  original  design  concept  for  such  a  log. 
shown  iri  Fig.  2.1  (Chilowski.  1932).  already  employs  the  sophisticated  four-beam  Janus 
configuration  still  iri  use  today. 

In  the  modern  instrument  used  for  this  study,  an  acoustic  transducer  mounted 
to  the  ship's  hull  simultaneously  transmits  a  short  pulse  of  acoustic  energy  (ping)  of  well 
defined  frequency  f0  along  four  beams,  each  inclined  30°  from  the  ship’s  vertical  axis, 
pointing  foreward.  aft.  port  and  starboard.  As  each  pulse  travels  down  its  beam,  it 
"ensonifies"  successive  volumes  of  ocean  along  that  beam  (Fig.  2.2).  Acoustic  targets 
within  the  erisonified  volume  may  scatter  a  portion  of  the  incident  energy  .  and  a  portion 
of  this  scattered  energy  will  be  directed  back  toward  the  transducer.  The  motion  of 
those  targets  relative  to  the  transducer  will  induce  a  Doppler  shift  Af  -  fR-f0  in  the  fre¬ 
quency  f|(  of  the  backseat tered  energy  received  at  the  transducer.  If  the  relative  velocity 
between  target  and  transducer  is  V,  and  the  speed  of  sound  is  c.  then,  to  first  order  in 

V,  c 

Af  —  (V.T-)  (2.1) 

< 

where  r  a  unit  vector  along  the  ilh  beam  To  this  order  of  approximation,  the  Doppler 
shift  is  therefore  a  measure  of  that  component  of  the  target  velocity  V,  which  lies  along 
the  acoustic  beam  direction 


Since  the  measurements  for  this  study  come  from  instrumentation  which  is 
somewhat  novel,  we  first  address  some  questions  regarding  the  accuracy  of  the  tech¬ 
nique.  The  importance  of  such  questions  may  be  readily  appreciated,  since  the  currents 
are  inferred  as  a  difference  between  two  directly  measured  quantities  dx0/dt  and  V,  each 
of  which  may  be  much  larger  than  the  current.  For  example,  to  measure  currents  to  an 
accuracy  of  1  cm  sec  from  a  ship  travelling  at  10  knots  (approximately  500  cm/sec), 
both  dx0  dt  and  V  must  be  measured  to  an  accuracy  of  0.2%.  After  a  brief  introduction 
to  the  shipboard  Doppler  system  in  Chapter  2,  the  Doppler  and  collateral  measurements 
are  examined,  questions  of  accuracy  are  explored,  and  considerations  required  in 
transforming  the  shipboard  referenced  measurements  into  an  Earth-fixed  reference  frame 
are  discussed  in  Chapter  3.  In  Chapter  4  Doppler  measurements  are  combined  with 
navigation  data  to  extract  the  currents.  Chapter  5  presents  a  comparison  of  measured 
currents  from  the  shipboard  Doppler  system  with  those  made  from  moored  current 
meters.  With  this  background  we  then  look  at  the  ocean  measurements  collected  during 
the  CODE  experiment  in  Chapters  6  and  7,  where  the  quasi-synoptic  field  is  mapped 
and  the  spatial  variability  in  the  mean  field  and  several  interesting  event  scale  features 
are  described.  We  find  that  the  synoptic  view  of  ocean  currents  shows  an  astonishingly 
rich  field  of  energetic  eddies  and  jets  which  actively  transport  upwelled  coastal  water 
offshore.  These  features  of  the  synoptic  current  field  occur  over  a  wide  range  of  scales. 
The  mean  field,  on  the  other  hand,  resembles  the  rather  smooth  earlier  findings  shown  in 
Fig.  1.1. 


circulation  (from  Huyer  1983) 


instrument  can  sample  different  locations  becomes  the  new.  less  restrictive,  constraint  on 
the  set  of  (x.y.z.t)  points  which  can  be  studied. 

Over  the  past  two  decades,  moored  current  meters  have  been  used  to  intensively 
study  the  current  field  iri  coastal  upwelling  regions.  Fig.  1.1  shows  a  schematic  picture 
of  tin  upwelling  circulation  which  has  been  derived  from  such  measurements  off  Oregon 
(Huyer.  1983).  Equatorward  wind  forcing  causes  offshore  transport  of  warm  surface 
water  in  a  surface  Ekman  boundary  layer.  To  replace  it.  cold  dense  water  at  depth  flows 
toward  the  coast.  This  cold  water  surfaces  in  a  band  near  the  coast.  The  alongshore 
currents  are  described  by  a  surface  intensified  equatorward  jet.  whose  vertical  shear  is 
geost rophically  balanced  by  cross-shore  density  gradients  induced  by  the  upwelling. 
Flowing  counter  to  the  direction  of  wind  forcing,  a  poleward  undercurrent  is  seen  at 
depth.  Models  of  this  circulation  usually  assume  that  gradients  in  velocity  are  small 
enough  that  advective  terms  can  be  neglected  in  the  momentum  balance,  and  that 
alongshore  variations  of  the  currents  are  much  smaller  than  cross-shore  variations 
(Allen.  1980). 

When  spatial  variability  can  be  determined  synoptically .  this  simple  picture 
derived  from  moored  instruments  gives  way  to  a  much  more  complex  one.  Satellite 
image'-  of  sea  surface  temperature  show  cold  upwelled  water,  rather  than  occurring  in  a 
simple  band  along  the  coast,  forms  very  complex  patterns  with  strong  alongshore  varia¬ 
bility.  I  pwellmg  centers,  patches  of  coastline  with  intense  upwelling.  have  been 
observed  in  many  locations  (see  Brink.  1983).  Tongues  of  cold  water  extending  several 
hundred  kilometers  out  to  sea  have  been  recognized  along  the  west  coast  of  the  T.S. 
(Bernstein,  el  al.  1972.  Breaker  and  Gilliand  1981.  Traganza.  el  al.  1981.  Kelly  1983). 
Because  these  features  are  hard  to  study  with  moored  instruments,  little  is  directly 
known  about  the  strength  or  structure  of  the  circulation  which  forms  them.  If  the 
features  are  associated  with  strong  currents,  they  may  represent  important  pathways  of 
exchange  between  coastal  and  offshore  waters. 

During  the  ( 'oast al  Ocean  Dynamics  Experiment  (CODE),  a  broad  range  of 
techniques  were  applied  to  the  measurement  of  ocean  variability  in  the  coastal  and 
offshore  waters  of  Northern  California  during  the  spring  and  summer  upwelling  seasons 
of  19N1  and  1982  (CODE  Croup  1983).  In  addition  to  extensive  moored  instrumenta¬ 
tion.  a  major  effort  was  made  to  resolve  the  spatial  structure  of  the  upwelling  Fields 
through  the  use  of  satellite  infra-red  (1R)  imaging  for  sea  surface  temperature,  aircraft 
mapping  of  sea  surface  temperature  and  meteorological  data,  a  large  program  of  drifter 
measurements,  as  well  as  shipboard  DAL  current  profiling.  This  study  reports  the 
results  of  the  DAL  program. 
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INTRODUCTION 


The  shipboard  Doppler  Acoustic  Log  (DAL)  to  be  discussed  in  this  study  is  a 
new  tool  for  the  measurement  of  ocean  currents.  By  providing  vertical  profiles  of 
currents  over  the  upper  150m  of  ocean  at  points  along  the  ship’s  path,  it  allows  tin- 
ocean  to  be  sampled  in  a  way  which  is  fundamentally  different  from  moored  current 
meters  or  drifters,  the  instruments  most  commonly  used  for  current  measurement.  This 
instrument  is  used  to  examine  the  circulation  in  a  region  of  strong  coastal  upwelling  off 
Northern  California,  where  it  reveals  an  unexpectedly  energetic  and  spatially  complex 
synoptic  flow  field. 

Consider  the  general  problem  of  measuring  the  current  at  the  water  parcel 
whose  coordinates  are  xw(t)  using  an  instrument  whose  location  is  x0(t).  explicitly  a 
function  of  time.  By  a  simple  identity 

dxw  dxCi  d(xw-  x(l) 

dt  dt  "  dt 


In  the  notation  which  will  be  used  throughout  this  study.  Eq  (1.1)  may  be  rewritten  as 


u(x0-r) 

where 


dx0 

dt 


-  V(r) 


(1.2) 


r  =  xw  -  x0 

u(xc-r)  =  dxw/dt 

V(r)  =  d(xw-x0)/dt 

The  current  u  at  the  measurement  point  x0-r  must  be  determined  from  the  sum  of  two 
terms.  One.  the  relative  velocity  V(r)  .  is  the  velocity  of  the  water  parcel  relative  to  the 
instrument  position.  The  other,  dxo/dt.  is  the  velocity  of  the  instrument  itself  with 
respect  to  the  Earth.  Moored  instruments  and  drifters  each  measure  only  one  of  these 
terms,  and  are  engineered  to  make  the  unmeasured  term.  dx0  dt  for  moored  instruments 
and  V  for  drifters,  negligible.  This  design  constraint  places  limits  upon  the  types  of 
variability  which  each  instrument  can  sample.  Although  ocean  currents  vary  in  all  three 
spatial  dimensions  as  well  as  in  time,  a  moored  instrument  samples  only  the  temporal 
variability  at  a  single  location,  while  a  drifter  samples  the  time  and  space  variability 
only  along  the  path  of  a  single  water  parrel.  For  an  instrument  such  as  the  DAL.  which 
directly  measures  both  V  and  dx0/dt.  these  constraints  on  possible  sampling  trajectories 
are  removed,  allowing  mapping  of  the  current  field.  The  speed  with  which  the 
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Table  3.1  :  Nominal  D\L  operating  parameters 

Table  3.2  :  Range  bin  geomem 

Table  3.3  :  Gyroscopic  measurements  of  the  ship’s  attitude  were  used  to 
correct  the  shipboard  profiles  of  relative  velocity  V(z)  for  the  effects  of  pitch  and 
roll  during  CODE  1.  This  table  shows  the  errors  (in  cm  sec)  which  would  have 
been  incurred  in  100  ping  averages  of  V ( z )  in  the  absence  of  such  corrections. 
Less  than  1  cm  sec  of  bias  is  introduced  into  the  horizontal  components  of  V(z) 
when  pitch  and  roll  compensation  is  not  performed.  The  added  noise  in  the 
estimates  is  likewise  small.  These  calculations  were  performed  on  data  subsets 
described  in  section  3.0  of  the  text. 

Table  3.4  :  Heading  Errors.  This  table  lists  several  common  errors  which  may 
be  present  in  gyrocompass  data.  Magnitudes  assume  an  operating  latitude  of  38° 
N.  a  course  due  north  at  10  knots,  and  a  maneuver  in  which  the  ship  reverses 
course  over  a  time  which  is  short  compared  with  the  Schuler  period.  With  the 
exception  of  the  rolling  error,  these  mat  be  corrected  either  at  the  gyrocompass 
or  in  post-processing. 

Table  3.5  :  Calibration  error  :i  and  misalignment  angle  66^,  calculated  from 
changes  in  apparent  current  which  accompany  changes  in  ship  velocity  using 
Eq  3.12.  The  large  value  of  <*>0,  during  CODE  1.  Leg  4  decreases  substantially 
for  subsequent  cruises:  this  drop  coincides  with  repair  and  realignment  work 
performed  on  the  gyrocompass  following  Leg  4. 

Table  5.1  :  CODE  2  moored  array.  S  (SS)  denotes  a  current  meter  string 
which  is  surface  (sub-surface)  moored.  \  ACM  (VMCM)  denotes  a  Vector 
Averaging  (Measuring)  Current  Meter. 

Table  5.2  :  Current  meter  DAL  comparison  statistics.  Whenever  the  DAL 
was  within  1  km  of  a  mooring,  comparison  of  DAL  and  moored  current  meter 
measurements  was  made.  Averages  and  variances  for  each  instrument  are  shown, 
as  are  covariance,  correlation  and  standard  deviation  of  the  difference  between 
measurements  from  each  instrument.  Also  shown  is  a  comparison  between  the 
only  closely  spaced  pair  of  moored  instruments  in  CODE,  a  VACM  and  a 
VMCM  both  at  10m  near  C3. 
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locations.  The  receiver's  internal  software  convert-  position  from  tii<  l.OR.W  grid  to  A 
and  L.  Tlie  system’s  advantages  lie  in  its  wide  (overage  and  continuous  availability. 

LORAN-C  fixes  are  subject  to  a  number  of  first  order  errors,  including  dual 
solutions  and  cycle  errors  (lane  jumps).  Dual  solutions  arise  when  only  two  time  delay 
measurements  (three  transmitters)  are  available.  Since  the  line  of  positions  which  yield 
a  given  time  delay  is  a  hyperbola,  and  hyperbolae  form  closed  curves  on  a  sphere,  two 
lines  of  position  will  intersect  at  two  distinct  points  in  general,  and  the  intersection 
corresponding  to  the  receiver  location  must  be  selected.  The  solutions  are-  often  well 
separated,  in  which  case  selection  of  the  proper  one  is  easy.  Cycle  errors  are-  discrete1 
offsets  in  the  time  delay  measurement  by  multiples  of  10  //sec.  They  occur  most  com¬ 
monly  when  the  signal  is  weak,  arid  are  caused  by  the  misidentificatiori  of  the  arrival 
time  of  the  pulse,  defined  by  the  third  rising  zero  crossing  of  the  100  kHz  carrier  within 
the  pulse.  If  the  wrong  zero  crossing  is  tracked,  the  time  delay  measurement  slips  by  an 
integral  number  of  carrier  periods.  Such  a  jump  usually  results  in  a  position  displace¬ 
ment  of  several  miles,  and  so  is  easily  delected.  Since  the  LORAN-C  time  delays  were 
recorded  along  with  the  Northstar  A  and  L.  cycle  errors  were  corrected  when  they 
occurred  by  correcting  the  time  delay  and  recalculating  the  associated  position  change- 
using  the  algorithm  of  Campbell  (19(18).  This  procedure  was  also  used  to  correct  for 
times  when  the  secondary  solution  was  erroneously  recorded 

More  subtle  errors  also  rnay  be  present  in  LORAN-C  position  data.  For  exam¬ 
ple.  the  time  a  pulse  takes  to  reach  the  receiver  depends  not  only  on  the  path  length  but 
also  ori  radio  wave  propagation  speed,  which  in  turn  depends  on  conductivity  along  the 
path.  Day  night  (Dean  1978)  and  land  sea  path  (Johler  et  al.  1956)  differences  may 
thus  be  important.  Another  concern  is  that  the  data  are  internally  filtered  by  a 
proprietary,  arid  thus  unknown,  scheme  before  being  output  by  the  receiver.  Any  lag 
produced  iri  position  during  accelerations  could  seriously  affect  the  usefulness  of  the  data 
for  determining  ocean  currents. 

The  most  direct  measure  of  accuracy  for  LORAN-C  fixes  comes  from  the  varia¬ 
bility  in  H.xes  taken  at  a  single  location.  Such  data  were  obtained  over  an  eight  hour 
period  while  the  ship  was  docked  at  \  erba  Buena  Island  between  cruises.  (An  interrup¬ 
tion  of  ships  power  aborted  the  data  collection  much  earlier  than  planned.)  The  rms 
noise  during  this  test  was  13  rri  in  latitude  and  27  m  in  longitude.  The  principal  axis  of 
variability  was  along  2-19'T:  the  major  and  minor  axis  rms  errors  were  29  in  and  8  m 
respectively.  No  significant  serial  correlation  was  found  between  successive  position 
errors  over  the  sample*  of  62  fixes,  spaced  8.5  minutes  apart. 


\  hese  accuracy  estimates  arc  likely  to  be  optimistic  for  our  purposes.  since  they 
do  not  include  contributions  from  receiver  motion  and  probably  undersample  long  term 
variability  .  Tin  question  of  position  lag  during  ship  accelerations  due  to  internal  filter¬ 
ing  was  examined  by  computing  lagged  correlations  between  the  ship  velocity  as  meas¬ 
ured  by  tin  acoustic  system  and  ship  velocity  inferred  from  LOHAN-C  fixes.  Any 
filter-induced  lag  would  show  up  as  ari  asymmetry  between  positive  and  negative  lagged 
correlations.  No  such  asymmetry  was  found  for  lags  separated  by  the  sampling  period  of 
70  sec.  Therefore,  for  our  analysis  we  shall  assume  that  the  LORAN  yields  position 
measurements  which  contain  errors  of  order  10  to  30m  and  that  these  errors  can  be 
modelled  as  white  noise,  i.e.  uncorrelated  in  time. 

3.3.  Thermistor 

A  manufacturer  supplied  thermistor  provided  ocean  temperature  data  at  the 
acoustic  transducer.  This  data  was  used  in  the  estimation  of  sound  speed  at  the  trans¬ 
ducer.  ('lay  and  Medwin  (1977)  write  for  the  speed  of  sound  (m  sec) 

c  14-19.2  -  4.6T  0.055T-  -  .00029T''  -  (1.34  -  0.010T)(S  35)  -  .016z 

where  T.  S  and  z  are  the  temperature  in  °C.  salinity  in  ppt.  and  depth  in  m.  Thus  a  rise 
in  temperature  of  i  "C  increases  the  sound  speed  0.3(/<:  an  increase'  in  salinity  of  1  ppt 
results  in  a  O.lT  increase  in  c  From  the  range  of  variability  seen  in  surface  temperature 
and  salinity  maps  in  the  ('OI)E  region  (  e.g.  Fleischbein.  Gilbert  and  Huyer.  1982) 
corrections  to  the  sound  speed  due  to  surface  salinity  changes  are  at  or  below  the  0.1(7 
level,  hence  negligible,  while  the  corrections  due  to  temperature  variability  can  be  1(7  or 
more,  hence  essential.  Snell's  law  implies  that  vertical  variability  in  c.  while  larger  than 
horizontal  variability,  does  not  affect  the  right  hand  side  of  (2.1)  for  the  horizontal  com¬ 
ponents  of  V  w  hich  are  of  i.  terest  . 

The  thermistor  data  also  provides  an  underway  map  of  near  surface  tempera¬ 
tures  concurrent  with  the  acoustic  data.  The  data  quality  was  estimated  by  comparing 
the  measurements  with  CTD  temperature  profiles  for  Leg  4  of  CODE  1.  Only  stations 
for  which  the  upper  water  column  was  well  mixed  were  used  in  making  the  comparison. 
The  thermistor  temperatures  showed  a  mean  error  of  0.4'C.  The  mis  variations  about 
this  mean  difference  decrease  with  time  on  station,  implying  that  the  thermistor  is  insu¬ 
lated  by  the  acoustic  transducer  head.  The  time  history  of  some  of  the  larger  differences 
indicates  an  exponential  decay  time  for  equilibration  of  order  20-30  minutes. 
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3.4.  Data  Valid  Flag/  Depth  of  Acoustic  Measurements 

To  exclude  grossly  erroneous  data,  acoustic  velocity  estimates  at  any  depth  were 
rejected  utdess  all  4  beams  showed  the  "data  valid"  flag  set  at  that  depth.  This  flag,  set 
In  the  Ametek  hardware,  is  an  indication  that  the  backscattered  energy  received  at  the 
transducer  was  sufficient  to  resolve  the  Doppler  shifted  peak  of  the  spectrum. 

fig.  3.1  shows  how  the  average  data  quality  over  a  12  hour  period,  as  indicated 
In  this  flag,  varied  as  a  function  of  depth  for  2  subsets  of  the  data  taken  in  water  at 
least  500  m  deep.  The  details  between  closely  spaced  vertical  profiles  can  vary,  but  in 
general  there  is  first  a  region  near  the  surface  for  which  the  data  quality  is  poor,  next  a 
broad  region  where  data  quality  is  excellent,  and  finally  a  region  in  which  the  data  qual¬ 
ity  diminishes  with  depth. 

Since  the  intensity  of  scattered  energy  falls  rapidly  with  distance  from  the 
scatterer.  it  is  not  surprising  that  data  quality  diminishes  with  depth.  It  is  conjectured 
that  the  near  surface  region  of  low  validity  is  associated  with  spectral  spreading  caused 
by  scattering  from  bubbles  or  directly  from  the  ships  hull.  A  remarkable  example  of 
how  scattering  strength  can  vary  in  unforeseen  ways  is  reported  by  Cochran  and 
Sameolo  ( 1983). 

3.5.  Bottom  Reflection 

Since  CODE  was  a  coastal  experiment,  acoustic  data  was  often  collected  in 
waters  for  which  the  bottom  depth  was  less  than  the  acoustic  range.  In  such  regions,  it 
was  necessary  to  exclude  signals  reflected  from  the  seafloor.  The  profile  of  received 
power  may  provide  a  straightforward  means  of  accomplishing  this,  but  such  measure¬ 
ments  proved  unreliable  during  the  experiment.  Screening  was  thus  performed  on  each 
average  profile  based  on  the  bottom  depth  estimated  from  the  ship's  position. 

Iri  the  CODE  region,  ari  efficient  computer  routine  was  developed  for  interpolat¬ 
ing  a  digitized  bathymetry,  using  LORAN-C  fixes  for  position.  The  accuracy  of  the 
predictions,  by  comparison  with  depth  sounder  (PTR)  measurements,  was  markedly 
improved  by  adjusting  measured  LORAN-C  positions  0.3'  north  and  0.25  east.  This  is 
consistent  with  a  comparison  of  satellite  vs.  LORAN-C  derived  positions.  After  adjust¬ 
ment  of  LORAN-C  positions,  the  rms  accuracy  of  the  depth  predictions  within  the  DAL 
depth  range  was  2.5  in.  or  less  than  1  range  bin.  Outside  the  CODE  area,  minimum 
depth  sounder  readings  over  a  time  interval,  generally  5  minutes,  were  read  manually 
and  keypunched  for  interpolation.  In  the  absence  of  depth  sounder  records.  LORAN-C 
positions  corrected  by  local  estimates  of  the  LORAN-C  offset  from  satellite  fixes  were 
plotted  on  depth  charts  and  depths  keypunched  for  interpolation. 


depth  (m) 


Fig  3.1  :  Depth  range  of  shipboard  measurements  is  shown  by  the  percent  of  pings  for 
which  the  profiler  electronics  indicated  valid  returns  at  each  range  bin.  Two  12  hour 
subsets  of  data  are  shown,  one  each  from  1981  and  1982.  Details  of  such  validity  profiles 
vary  with  many  factors  including  location  and  time  of  day. 


Si  lie  t  l)ul  i  on  i  reflect  ion  is  expected  to  be  much  stronger  than  ocean  backscatter. 
reflections  from  transducer  sidelobes  will  degrade  the  signal.  To  reject  all  possible 
sidelobe  reflections  from  the  bottom,  the  profile  was  regarded  as  valid  only  for  range  bins 
whose  maximum  range  is  less  than  the  water  depth.  D.  minus  t fit-  transducer  depth.  z<„ 
From  Table  3 . .  relating  range  biri  geometry  and  depth,  the  relative  velocity  profile  was 
regarded  as  \alid  at  depth  z  only  if 
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or  for  the  nominal  values  of  the  parameters 


/  0.S7I)  9  meters. 


Because  of  tin  need  to  eliminate  bottom  reflection,  near  bottom  currents  cannot  be  pro¬ 
filed. 


3.6.  Acoustic  Profiles  of  Relative  Velocity  in  Ships  Coordinates 

How  accurate  art  acoustic  measurements  of  relative  velocity  V  made  from  a 
moving  ship?  They  will  contain  errors  from  a  wide  variety  of  sources,  arid  we  may 
expect  large  variability  between  individual  profiles.  Variability  arises  from  measurement 
noise  due  to  finite  bandwidth  in  the  transmitted  acoustic  pulse,  as  well  as  finite  time  for 
frequency  resolution  within  a  range  bin.  Variability  occurs  in  the  ocean  at  scales  below 
our  resolution,  from  small  scale  turbulence  within  a  range  bin  or  horizontal  shears  at 
scales  smaller  than  the  separation  between  acoustic  beams.  Even  if  the  scales  of  ocean 
velocity  were  fully  resolved,  additional  variability  is  introduced  by  limited  sampling  of 
non-passive  motion  (swimming)  of  the  objects  scattering  acoustic  energy  and  by  multiple 
scattering  events.  Accelerations  of  the  ship,  in  response  to  ocean  waves  or  under  inten¬ 
tional  control  from  t  fie-  bridge,  also  add  variability  to  the  relative  velocity  . 

fig.  3.2  shows  histograms  of  the  variability  in  measurements  of  F(z)  and  P(z). 
the  foreward  and  port  ward  components  of  V.  at  constant  z. the  depth  in  ships  coordi¬ 
nates  (Chap  2).  The  data  come  from  1-minute  segments  (100  pings)  over  which  the 

ship  speed  was  being  held  steady.  Despite  the  fact  that  the  data  are  taken  over  a  short 
time  interval,  the  mis  variability  in  individual  estimates  F  and  P  at  fixed  z’  were 

rrns  variability  on  station  at  10  knots 

in  f  38.9  cm  sec  36.3  cm  sec 

in  P  23.0  cm  sec  17.3  cm  sec 

The  size  of  this  variability  is  nearly  independent  of  the  range  bin  considered. 
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Fig  3.2  Variability  in  acoustic  measurements  of  relative  velocity,  shown 
histograms  of  the  deviations  of  single  ping  measurements  from  100  ping  average. 


Autospectra  of  this  high  frequency  variability  arc  shown  in  1  ig.  3 . 3  for  3  sets  of 
conditions.  In  eaefi  case  tiie  ship  is  steaming  at  "constant "  speed  under  control  from  the 
bridge.  The  spectra  show  that  the  variability  in  acoustic  estimates  over  periods  up  to  1 
minute  is  composed  of  a  white  noise  background  plus  a  broad  spectral  peak  at  periods  of 
order  3  to  10  seconds.  The  frequency  of  the  spectral  peak  is  higher  when  the  ship  steams 
into  the  waves  and  lower  when  steaming  away  from  waves.  These  characteristics  suggest 
that  the  peak  is  due  to  accelerations  of  the  ship  by  the  surface  wave  field.  This  is  con¬ 
firmed  by  the  high  coherence  found  iri  the  3  to  10  second  bands  of  the  cross  spectrum  of 
F ( 7  )  and  P(z')  with  the  pitch  and  roll  angle  of  the  ship  (Fig.  3.1).  The  contributions 
to  the  rms  variability  from  measurement  of  wave-induced  ship  accelerations  and  from 
the  white  noise  background  error  level  are 
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Clearly  the  majority  of  variability  at  periods  up  to  1  minute  arises  from  sensing  of  actual 
wave  induced  ship  motions  rather  than  from  white  noise  inaccuracies  iri  the  measure¬ 
ment  process  itself. 

Nonetheless,  before  they  can  be  used  to  infer  properties  of  the  ocean.  th<  meas¬ 
urements  must  be  filtered  to  reduce  both  the  white  noise  and  wave  induced  variability. 
The  success  that  any  proposed  filter  will  have  in  reducing  this  variance  can  be  calculated 
from  the  spectrum  of  the  noise  (Fig.  3.3)  and  the  transfer  function  of  the  filter.  Fig.  3.5 
shows  the  rms  noise  as  a  function  of  filter  length  N  for  a  block  averaging  filter.  The 
variance  is  reduced  rapidly  with  increasing  N  up  to  N  50.  then  more  slowly  for  higher 
\.  This  is  because  tin  variability  due  to  accelerations  in  the  wave  field  is  highly 
coherent  in  time  and  is  reduced  approximately  as  N  '.  The  white  noise  component,  while 
smaller,  decreases  only  as  N  1  *  and  so  persists  longer. 

We  see  then  that  by  regarding  the  individual  estimates  F  and  1’  as  time  series  at 
each  range  bin  and  averaging  over  many  pings,  the  high  frequency  variability  due  to 
waves  and  white  noise  inaccuracies  can  be  reduced.  Block  averaging  over  N-  100  sam¬ 
ples  leads  to  averaged  estimates  contaminated  by  1-2  cm  sec  rms  noise  from  the  residual 
effects  of  white  noise  measurement  errors  and  wave  induced  ship  motion.  As  N  is 
increased  further,  the  noise  from  these  sources  decreases  as  N  1  *. 
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Fig  3.3  :  Autospectra  of  ping-to-ping  relative  velocity  measurements  for  3  subsets  of 
data  taken  at  nearly  constant  shipspeed  while  (a)  steaming  into  seas,  (b)  on  station,  and 
(c)  steaming  with  following  seas.  Solid  line  is  foreward  velocity  F,  dashed  line  is 
portward  velocity  P.  Note  peak  at  surface  wave  encounter  frequency  rising  above  a 
white  noise  background. 


Fig  3.4  :  Coherence  of  F  (solid)  and  P  (dashed)  with  pitch  ^  and  roll  p  of  the  ship.  Data  in  left  panels  taken  while  on 
station  (same  as  middle  panel,  Fig  3.3);  data  in  right  panels  taken  while  steaming  into  seas  (  same  as  upper  panel, 
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Fig.  3.5  Expected  residual  wave  induced  and  white  noise  in  acoustic  measurements 
of  relative  velocity  after  block  averaging,  as  a  function  of  averaging  length. 


3.7.  Rotation  to  Geographical  Coordinates 


To  properly  vector  average  the  relative  velocity  measurements.  they  must  first 
he  transformed  from  t  h<-  ship's  coordinate  system  (x  .y  .7  )  to  geographical  coordinates. 
Let  us  define  a  (quasi)  geographical  coordinate  system  with  origin  at  the  ship's  trans¬ 
ducer  and  coordinate  axes  x.  y  and  z  pointing  east,  north  and  up  respectively.  We  shall 
describe  the  transformation  between  geographical  and  ships  coordinates  by  a  series  of 
three  rotations  involving  the  angles 

(I.  the  ship's  heading,  measured  from  y  (North)  to  x"  (Fore) 

p.  the  roll  angle,  positive  for  port  side  elevated 

<:• .  the  pit c  1 1  angle,  positive  for  bow  elevated 
measured  relative  to  rolled  coordinates 

Any  arbitrary  vector  whose  representation  in  shift  coordinates 
representation  in  geographic  coordinates  given  by 

A  RJR  ,RCA 


is  A  will  have  a 
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Thus  if  tin  relative  velocity  vector  measured  in  shift  coordinates  is  V  .  the  firm  aver¬ 
aged  relative  velocity  in  geographic  coordinates  will  be 


V  R„R  R  V 


Fq  (.'!..''>)  emphasize*  the  fart  that,  because  the  transformation  between  the  two  systems 
is  time  dependent,  we  must  transform  the  data  to  geographic  coordinates  before  vector 
averaging  to  obtain  correct  averages.  Otherwise,  correlations  between  the  shifts  orienta¬ 
tion  and  its  velocity  will  bias  the  resulting  average  and.  in  turn,  the  currents  calculated 
from  those  averages.  Since  the  transformation  is  performed  using  measured  values  of  0. 
P  and  C.  the  errors  in  measuring  these  angles  will  affect  inferred  average  velocities.  Tin- 
analysis  of  the  next  two  sections  suggests  (hat  the  effects  of  pitch  and  roll  are  not 


important  for  tin-  ('OI)L  data  set.  tint  that  error'  i i,  ships  heading  measurement s  might 
be. 


3.7.1.  Pitch  and  Roll 


\<  com. ting  ior  tin  eflect  of  pit 1 1 1  and  roll  of  tin  slop  on  measured  profiles  of  the 
relative  voiion  vector  reimire-  c  oii'idi-rai  ton  of  two  effects..  In  tin  first  place,  t  he  com- 
pomi.t-.  ci'  \  riiii't  be  rotated  to  level  coordinates.  In  addition,  the  location  of  the 
tin  d-i;ri  !ii‘  /  .  1 1 1 1 ; '  ’  a  l  sc '  t  n  t  rati' I  or  med  .  Thus 
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dP  dz 
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sinpdH  dz 

dH  dz 

taii'ciufeospsiriodF  dz-sinpdP  dz) 

i-  tin  term  which  must  be  added  to  the  profile  in  ships  coordinates  to  correct  for  dis¬ 
placement  of  the  measurement  depths  from  their  nominal  positions.  It  is  calculated  from 
tin  bean:  equations  2.7  and  2>  tinder  the  assumption  that  horizontal  shears  are  negligi¬ 
ble  (ompared  to  vertical  shears  The  full  error  vector  E  is  the  velocitv  overestimate  if 
the  profili  I-  riot  corrected  for  pitch  and  roll. 

E  ll  R,  R,  iV  iz  )  R,  Rf  i'V 


I  ( 1  cos o  )  -  Hstnc 
I’ll  cosp)  ■  1  sirtpeose  •  Hsinpcoso 
11(1  cosocosp)  fsinocosp  Psinp 


R(R^v 


(3.6) 


I  riles'  pitch  and  roll  compensation  i'  performed  separatelv  for  each  profile,  variability  in 
E  will  contribute  noise  to  estimates  of  V(z)  as  the  ship's  attitude  changes  from  profile  to 
profile.  More  important  Iv.  if  E  has  a  non-zero  mean  value  over  the  averaging  period 
used  this  mean  value  will  bias  the  estimate  of  V(z). 


"V - "W - T”  '  T 


21) 


Table  3.3  :  Effect  of  Pitch  o  and  Roll  p  angles 
on  Relative  Velocity  Components  (cm,  sec) 


Depth  =  30  m 


Meaming 


On  Station 


mean 

std  dev 

max 

min 

mean 

std  dev 

max 

min 

Measurements 

F 

500. 

1 1.2 

523.2 

480.1 

19.8 

31.8 

92.5 

-28.7 

F’ 

2.4 

3.1 

1  1.8 

-1 1.8 

1.1 

19.0 

41.1 

-37.0 

H 

-1.6 

1.6 

3.1 

-5.9 

19.2 

2.1 

26.6 

15.3 

Depth  change  error 

Tf 

.0 

0 

.0 

.0 

.0 

.0 

.0 

.0 

CP 

(1 

.0 

.0 

.0 

.0 

.1 

.0 

.0 

cH 

.0 

.0 

.0 

.0 

.0 

.0 

.0  : 

.0 

Kertificat  ion  Error 

F  ore 

F |  1  - coso  i 

1 

.0 

.3 

.1 

.0 

.0 

.0 

.0 

Hsino 

- 

‘» 

1.4 

.3 

.6 

.2 

1.0 

2 

Port 

P(  1  -  cos p  I 

.0 

.0 

.() 

.0 

.0 

.0 

2 

-.i 

Fsinosinzi 

-i 

.0 

.0 

.2 

.0 

.0 

.0 

.0 

Hcososinp 

-.i 

.1 

2 

-.4 

.0 

.3 

1.0 

-1.2 

Heave 

HI  1  -cosocosp ) 

0 

.0 

.0 

.0 

.0 

.1 

.0 

Fsinocosp 

■4.1 

1  .2 

:.i 

.5 

-.2 

.3 

.2 

-1.4 

Psinp 

.1 

1 

.3 

.0 

.6 

.6 

3.0 

.0 

Set  error  icm  scci 

in  F 

.  * 

2 

1.6 

.4 

.6 

.2 

1.0 

.2 

in  P 

-.1 

,i 

2 

-.5 

.0 

.3 

1.0 

-1.1 

in  H 

4.2 

1.2 

7  2 

.6 

.4 

.  < 

3.1 

-.9 

Table  3.3  :  <■  v  rosropir  measurements  of  the  ship's  attitude  were  used  to  correct  the  ship¬ 
board  profiles  of  relative  veiocitv  Viz)  for  the  effects  of  pitch  and  roll  during  CODE  1.  This 
table  shows  the  errors  (m  cm  sect  which  would  have  been  incurred  in  100  ping  averages  of 
V ( z )  in  the  absence  of  such  corrections.  Less  than  1  cm  sec  of  bias  is  introduced  if  to  the 
horizontal  components  of  Viz)  when  pitch  and  roll  compensation  is  not  performed.  The 
added  noise  in  the  estimates  is  likewise  small.  These  calculations  were  performed  on  data 
subsrrs  described  in  section  3.6  of  the  text 


During  ('ODE  I  the  pitch  ami  mil  of  the  ship  were  measured  gyroscopically  dur¬ 
ing  all  acoustic  profiling  as  described  in  Section  3.1  These  data  were  used  to  calculate 
tin-  terms  in  E  which  are  introduced  when  pitch  and  roll  corrections  to  the  profile  of 
relative  velocity  are  neglected.  Results  art-  presented  in  fable  3.3  for  two  data  subsets 
during  which  t  hi  ship’s  heading  0  was  held  constant,  one  w  hile  the  ship  was  steaming, 
tin  other  while  the  ship  maintained  position.  Each  data  subset  consisted  of  (i(l  blocks  of 
100  pings  each  (approximately  1  hour)  and  correspond  to  the  data  whose  cross-spectra 
were  presented  in  fig.  a.-).  The  results  indicate  that,  while  the  long  term  average  of  ship 
vertical  velocity  shows  a  substantial  bias  of  order  3--1  cm  set  when  pitch  and  roll  com¬ 
pensation  is  not  performed,  the  long  term  bias  in  the  horizontal  velocities  is  less  than  1 
cm  sec.  Table  3.3  indicates  that  the  major  contribution  to  the  bias  in  the  fore  com¬ 
ponent  arises  from  rectification  of  the  heave  component  through  the  term  HsiriO.  while 
the  major  contribution  to  the  bias  in  vertical  velocity  comes  from  the  Fsinocosp  term. 
The  effect  of  range  bin  displacement  on  long  term  averages  was  in  general  very  small. 
The  rms  size  of  the  additional  noise  in  a  100  ping  average  profile  due  to  variability  in  E 
is  of  order  0.3  cm  sec.  which  is  small  compared  with  the  1  to  2  cm  sec  noise  from  high 
frequency  variability  which  remains  after  filtering,  as  discussed  in  the  previous  section 

These  results  indicate  that,  for  a  well-riding  vessel  sue!)  as  t fie  R  \  Wecoma. 
averages  of  acoustic  profiles  can  be  calculated  with  minimal  error  in  horizontal  velocity 
estimates  by  treating  the  data  as  though  they  were  measured  in  a  level  plane,  ignoring 
corrections  for  pitch  and  roll. 

3.7.2.  Heading 

The  third  reference  angle  which  must  be  determined  before  the  relative  velocity 
profile  cati  be  rotated  to  geographic  coordinates  is  the  ship’s  heading.  Surprisingly,  it  is 
measurement  of  this  angle,  not  pitch  and  roll,  which  introduces  the  largest  uncertainty 
into  the  data  required  for  determination  of  currents.  This  comes  about  because  a  small 
measurement  error  in  0.  the  heading  angle  between  No-  h  and  the  fore  direction,  can 
lead  to  large  errors  in  the  inferred  geographic  components  of  relative  velocity,  and  there¬ 
fore  to  large  spurious  currents.  If  the  measured  heading  is  0  0  -  <1>f/.  where  b 0  is  the 
measurement  error,  then  the  geographic  components  of  relative  velocity  computed  using 
Eq  (3.2)  will  be  in  error  by 

l  (/I  I  (/)  I  (/)(cosM  1)  ■  V(z)sinM 

\(zl  \(z)  I  (zlsmZif?  \(z)(cos  t>0  1) 

for  a  '-hip  steaming  north  at  V  .’>00  cm  sec  (  10  knots.  ).  a  heading  error  frV  of  1 
would  v  ield  errors  of  (*.7.-0. 1  )  cm  sec  in  the  relative  velocities  (l  .Y).  and  so  add  errone- 
mi'-  contributions  of  the  same  size  to  the  currents.  Since  the  error  is  dominated  by 


am:  seek  l  Ik  values  of  n.  <1  which  product  I  lie  minimum  tola!  squared  mtsfil  between 
the  data  ami  tin  model 


V'ui  lit  d,  r. 


( -i .  :t  i 


Not.  that  <i  U  must  be  treated  as  a  measttrerneiit .  This  T  because  t  he  best  fit  will  riot, 
it.  general.  pa's  through  d  0  since  all  fixes,  including  x , contain  ineasureineiit  error. 
Minimizing  *'  with  respect  to  ti  and  d  gives 
\  N  N 

(\-l)£d;t;  Ef,  E'l, 

»■  - H - — — -  i-'-3*-' 

(\-  1 1'Ti  -  (  t  I 


V.lt  V; 
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This  least  squares  fit  tin*’  procedure  can  be  viewed  equallv  well  as  a  filtering  procedure 
This  nun  l ><■  seen  bv  rewriting  (4.32)  in  the  form 


yv.d 


(4.341 


(\1M  Du. 


|)E'/'  (S't 


(4.33) 


Tin  -  do  not  depend  on  the  data  values  d.  and  rnav  be  viewed  as  filler  weights.  If  the 
sampling  occurs  at  equal  interval' 


t  j  At 


(4.36) 


At  \(\  1  ) (  N  -2) 

di'tituting  |1.37|  and  (4.291  into  (4.34)  we  find 


1(1 


Nolr  dial  1  he-  n  1 1 e> r  is  suimicirK  .  w  w  'I  hi-  mean'  that  w (  need  not  have 
required  the  current  to  be  consiam  over  the-  filtering  interval.  I  he-  weaker  condition 

—  ( u  .  •  u  u  (4.25) 

2  -  -  J  ‘ 

suffices  for  (-t.l(i)  cl  *c<j.  to  apply. 

This  optimum  filter  is  a  significant  improvement  over  simple  averaging,  since 
the  rrtis  measurement  noise  in  the  estimated  current  decreases  as  \1  “  *.  compared  with 

tin  rate  \1  1  obtained  for  simple  averaging.  The  simple  average  has  a  clear  meaning  - 

rather  than  calculating  the  current  between  sequential  fixes,  average  \  over  (2\1  -  1)  fix 
intervals  and  use  only  the  endpoints  ji:.M  to  determine  tin-  current.  The 

optimum  estimate  also  lias  a  simple  physical  interpretation.  As  we  show  below,  it  is 

identical  to  making  a  least  squares  fit  to  the  slope  with  time  of  the  displacement  of  a 

water  parcel  at  the  measurement  depth 

4.5.  Least  Squares  Drift 

During  a  single  measurement  interval  At  a  water  parrel  at  the  measurement 
depth  will  drift  due  to  the  currents  by  an  amount 

11. At;  (X;  x,  ,)  \  .At  f,.  ( 4 . 20 ) 

VS  e  initially  assume,  as  before,  that  the  current  is  constant  over  the  water  parcels  sam¬ 
pled  during  the  estimation  interval.  Then  the  total  drift  from  time  t(1  to  time  tj 

d.  d.  x,  x,  -  £\,At  (4.27) 

;  i  i-  t 

will  imreasi  linearly  in  time  with  a  slope  given  by  the  current.  Thus  the  current  can  be 
estimated  by  using  measurements  (x,).{\  ,}  to  calculate  the  cumulative  drifts  and  then 
finding  tin  line  which  best  fits  them.  As  we  shall  see.  additional  useful  information  can 
bi  obtained  by  adopting  this  point  of  view. 

Let 

i  x  ■  x  |  \  At  (4.28) 

bi  tin  measured  drift  between  the  i  I'  and  i"  fixes.  Accumulate  this  in  a  total  drift 
from  the  start  of  the  estimation  to  the  j"  fix  in 

d  d  Vc  (4.29) 

i  i 

W  »  may  define  without  loss  of  generality  t,.  the  initial  time,  and  dfl.  the  measured  drift 
at  t  .  to  he  0.  Since  the  true  current  is  assumed  constant,  we  model  the  true  drift  as 


d  lit  •  d 


(4. HO) 


\-.-u  tump  that  the  efforts  of  Chapter  •’’>  wore  successful  and  the  bias  of  the  measurements 
An  (l.  I  in  bia-  in  I  Ik  filtered  estimate  i- 

M 

A  n  it  I  V  w  1  ].  (-1.17) 

\: 


Zero  bia-  i-  obtained  if  the  weight-  -J\\  l  [  satisfv 
M 

E  "  1-  (4. IS) 

M 

With  ( -1.  Hi )  and  (  1 . 1  s  1.  t  lie  mean  square  error  i- 
M  M 

An  "  E  E  w  ^ 11  .  Z* ti ,  _  (-1.19) 

M  ►  \1 

I  he  filter  weight-  art  determined  In  minimizing  this  mean  square  error,  subject  to  the 
constraint  H.lM  which  tan  be  enforced  with  t  Ik  aid  of  a  Lagrange  multiplier  A.  Letting 

M 

I  An'-  A(  E  w  •)  (4.2(1) 

J  M 

and  requiring  ell  dw,  (I  for  \1  n  \1  yields  the  2 VI  -  1  equations 
M 

E  w  •  A  u, .  ;Au. .  n  A  (4.21 ! 

M 


W  ith  the  constraint  (4. IS),  these  equations  define  at  most  one  solution  for  the  2M  -  1! 
variable-  (w  |  and  A.  Substituting  the  covariance  from  (4.12)  into  (4.21)  and  solving 
subject  to  the  constraint  (4. IS)  yields 
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_ a _ 

(  M  1 )  ( 2  M  -  I  )(2M  -  3) 


(\i  •  n1' 


M  j  '  M 


(4.22 


and 


/\ 


_ (> _ 

(  \l  -  1  K2.M  1 ) ( 2 M  -  .*{) 


which  i-  also  the  mean  square  error  of  the  estimator 


An"."  A  ft;  - — -  An' 
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(4.2 


(4.24 


as  mat  be  verified  from  (4.19). 


The-  rrn-  noise  in  u  .  ah  me  it \  id na!  t—1  mint e  < >f  i  h<  <  u rrent .  is  1  lien 


\  2r  , 

A; 


This  noise  is  epiite  large  :  tor  At  70  mi  .?  |.  s,  <  n  mc  and  10  cm  see  for  the  ross-sheire 
and  alongshore  i  urre-m-  re  -  pee  1 1\ e  ix 

''I  nee  ihi-  noise  o  proport  loual  lo  1  At  ai  iibxiou'  met  Mod  lor  improving  the 
aeenraev  «e)llle!  lx  III  Use  lunge'  lime  inlerval-  Ixtxxe-ci  fixe-  Increasing  At  to  \  At 
result',  m  an  \-|old  eh-e  re  as«  n  tin  rim  erro:  r  Hot  - 1 1 1  ■  pi  \  increasing  tin  time 
between  fixe"  Is  e-cim  x  a  n  I  '  In  axe  raging  \  <>«  j  ui  <  e  - 1 1 :  n  mia-'irei:  mer  the  shorter  Itilerxa 

At.  as  max  lx  seen  Iron  17  W  nn  . me  me  r-  are  a-a.lahn  m  Hie  -tioriiT  time 

interval  At.  sinip.t  averaging  thu-  e  <  .r  re -|>one’-  le.  throwing  awax  tin  n.iormai  iex,  pm 
xided  bx  intiTineeiiale  fixe-  .me:  using  emix  .  ii<  mitia.  and  final  fixe-  (  >i  the  other 
hand.  Ox  taking  tin  -latm'i'  o'  tin  Mgna:  am!  error'  m.o  aienunt  xx  e  ear  improxe  the 
fl  It  (ring  hex  one!  t  ne  sj  ii  ipii  ,ni  [aging  x1  a  n  e  m  1  (>"i  n  t  estimator 

4.4.  Optimum  Estimator 

liie  insight  gamed  from  spixira!  andlxsi-  ol  the  raxx  measurements  can  lx  used 
to  produee  an  estimator  ii  lor  l  In  eerremi  it  which  optimal!)  reduces  the-  noise  due  to 
errors  in  the  measurements  u  It  is  natural  to  const  rurt  this  estimate  from  a  linear  torn 
hination  of  the  \  AM  I  c  losest  measure  merits  centered  at  u,. 

M 

ii  V  xx  ii  (4.11) 

v 

Ihi  filter  weights  -|xx  .  are  lo  be  determined  l>>  minimizing  the  expected  mean  square 
error  in  ii  for  ,i  fixed  filler  lengtli  N  2M  -  I.  The  error  in  the  estimate-  (4.14)  will  be- 
v  M 

cn  M  xx  cu  M  xx  u  ni  .  (4.1-7) 

\i  \t 

The  tirsl  term  m  (  t  1  "> )  m  the-  err  or  due  to  noise  f'u  in  1  lie  measure  men  t  s.  xx  bile  the  term 
iri  brackets  is  the  error  which  would  arise  simplx  from  filtering  the  field,  even  iri  the 
absence  of  measurement  noise  (  learlx  it  is  not  possible'  to  c  hoose  a  single  set  of  weights 
which  minimizes  this  latte-r  term  for  arbitrarx  structure-  in  the  true  field  u,.  Since  the- 
analxsis  ol  the  preceding  section  suggested  that  the  geophxsical  signal  is  contained  in  the- 
low  frequcncx  portion  of  the  spectrum,  we  tentative!)  assume  that  u  is  constant  over  the 
filtering  interval  (xxe  shall  relax  ibis  condition  later).  Then 


•>  “ 


To  t  t*st  this  ititerpretai  ion  of  the  spectral  peaks,  we  consider  a  simple  model  for 
the  signal  and  noise  statistics.  The  error  spectrum  is  determined  by  the  correlation  func- 
l  ion  of  t  he  errors 

M 

it  M 

l.q  !  !.(>•  relate-  C 11  to  tin  errors  in  the  measured  quantities.  C\.  M  The  discussion  in 
('hapter  3  suggest-  that  it  is  plausible  to  model 

f  \  r  \  c  ~  t 

(  \  t\  r  \  t  (1.11) 

z  \  z  V  (I 

li  tin-  model  is  valid.  then  the  noise  in  the  current  measurements 

4^  (  (2  -  (-^— )*  )<\,,„  -  )  M-l-'! 

A  c , 

so  that  »i  should  observe  a  measured  spectrum. 

hr  '  „ 

> .  ( f )  S,„.  -  -  —  sin'g  fAt  -  2rr\"At  ( 4 . 1 '{ ) 

Hep  rest  tit  at  iv e  values  of  o  .'{()  m.  r  \  ~  2  cm  sec  and  Ai==  70  sec  are  cited  from  the 
result-  of  (’hapter  Errors  in  the  position  are  the  dominant  source  of  noise  over  most 

I  A  N  At 

of  tin  spectrum,  i.e  for  f  - sin  ( - 1  ■  10  Hz.  We  therefore  expect  that. 

'  At  2r  , 

il  tin  model  statistic-  are  correct,  the  spectrum  should  be  proportional  to  sin*7rfAt  over 
tht  noise  dominated  frequencies. 

Iii  f  ig.  -1,,''.  we  plot  S u u { f )  sin*7rfAt.  The  constancy  of  this  quantity  for 
f  1  2  •  10  Hz  supports  both  the  identification  rtf  the  high  frequency  peak  as  the 

noise  portion  of  the  spectrum,  and  also  verifies  the  statistical  model  (-1.1  1  )  over  this  most 
encrget  i<  part  of  the  inferred  noise.  Below  10  "  Hz.  u  contains  either  signal  from  the 
current  field  or  non-white  noise.  In  the  absence  of  other  information,  we  adopt  the 
former  possibilitv  as  a  working  hypothesis 

4.3.  Endpoint  Estimator 

Application  of  (1.13)  to  the  data  in  Fig.  4.3  allows  ari  estimate  of  the  mis  noise 
in  position  measurements 

re,  43  m  in  the  cross-shore  direction 

20  rrt  in  the  alongshore  direction 


jr  it  r  ri. 
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Fig.  4.1  sIiiihs  represent  at  i  v  «•  time  series  of  I  In-  measured  quantities  \  Ax,  At,,  and  ur 


4.2.  Spectrum  of  the  Measurements 

The  fx-to-fx  current  measurements  (u  }  are  a  time  series,  whose  spectrum  S  -(f) 
mat  tx  evaluated  via  a  discrete  Fourier  transform.  If  the  fixes  are  equal!)  spared  so 
t  hat  At  At  for  all  i. 
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(4.7) 


Then  from  (4.4).  assuming  signal  and  noise  are  uncorrelated . 


S,J  n  ■  suu(n  -  >wu»o  (4.9) 

i.e  .  t  he  spectrum  of  the  measurements  is  the  sum  of  the  noise  and  true  Field  variability 
spectra 

Figure  4.2  shows  the  spectrum  S  (f)  of  the  fix-to-fix  current  measurements  u 
for  a  subset  of  the  data  described  iri  Section  U.(>.  The  solid  curve  is  for  measurements  of 
the  current  component  along  SIT^T.  the  alongshore  direction  in  (’OI)E:  the  dashed  curve 
is  for  cross-shore  currents.  The  variance  in  the  measurements,  given  by  the  area  under 
the  spectrum,  is  larger  for  the  cross-shore  currents  than  for  the  alongshore  currents. 
Moth  components  show  similar  spectral  shapes,  characterized  by  a  sharp  peak  at  low  fre¬ 
quencies.  a  broad  rise  at  higher  frequencies,  and  a  spectral  gap  in  between.  It  is  tempt¬ 
ing  to  identify  the  low  frequency  peak  as  the  geophysical  signal  Suu ( f )  and  the  broad 
high  frequency  peak  as  the  contribution  of  the  noise  spectrum  S^Jf).  The  spectral  gap 
between  these  two  would  be  encouraging,  indicating  that  the  energetic  part  of  the  noise 
could  be  filtered  with  minimum  loss  of  the  geophysical  signal. 


Chapter  4 


ESTIMATION  OF  CURRENTS 


4.1.  Introduction 

\N  now  address  t  he  problem  of  obtaining  the  best  estimate  of  the  current  field 
using  the  data  described  in  Chapter  3.  In  Chapter  1  we  wrote 

u(x0-r)  -  -  V ( r )  (4-1) 

dt 

where  xk.  is  t he  location  of  the  instrument,  r  is  the  location  of  the  measurement  relative 
to  the  instrument,  and  u.  dx0  dt  and  V  are  the  current,  instrument  velocity  and  relative 
velocity  respectively.  Our  data  set  consists  of  measurements  of  consecutive  fixes  x0  and 
of  the  average  vertical  profile  of  relative  velocity  V(z)  between  fixes. 

First  consider  the  problem  using  measurements  from  a  single  depth,  with  each 
velocity  component  treated  sep  -ately.  Using  tildes  to  identify  measured  quantities  and 
hats  to  identify  estimates,  we  define 

t;  the  time  of  the  ith  position  fix  (4.2a) 

At,  -  l,  ■  I-,.. ,  (4.2b) 

X|  x0(t,).  the  true  position  at  time  t,  (4.2c) 

Xj  -  x,  -  .  the  measured  position  at  time  t;  (4. 2d) 

l, 

Yj  .£  —5—  f Y(zk)dt  .  the  true  average  V  between  fixes  (4.2e) 

A,i 

t , 

Vj  r  — —  f  V ( z k ) d t  -  V,  •  M’j  .  the  measured  average  Y  between  fixes  (4.2f) 
At  ,  J  t 

where  zk  is  the  depth  of  the  klt  range  bin  and  bxr  f>\  ,  are  the  errors  in  the  measure¬ 
ments. 

Each  pair  of  fixes,  together  with  the  average  relative  velocity  between  them, 
constitutes  a  measurement  of  the  average  current  at  zk 

u,  -  -^j-(Xj  Vi)  ^  A,  (4-3) 


f  rom  (4.1)  and  (4.2)  then 


Table  3.5  :  Calibration  and  Misalignment  Errors 


Samples 

Misalignment 
SO o  (degrees) 

Calibration 

3 

CODE  1.  Leg  4 

161 

-1.4 

-.006 

Leg  5 

153 

-0.2 

-.012 

Leg  7 

138 

0.1 

-.010 

CODE  2.  Leg  4-6 

174 

0.4 

-.009 

Leg  9 

174 

0.0 

-.010 

Table  3.5  :  Calibration  error  3  and  misalignment  angle  SO^  calculated  from  changes 
in  apparent  current  which  accompany  changes  in  ship  velocity  using  Eq  3.12.  The  large 
value  of  S80  during  CODE  1.  Leg  4  decreases  substantially  for  subsequent  cruises:  this 
drop  coincides  with  repair  and  realignment  work  performed  on  the  gyrocompass  follow¬ 
ing  Leg  4. 


Obtaining  current  est  i mat <•>  u  find  v  i-  discussed  in  (  hap.  1.  W  illi  thesi  est  i- 
mates  and  tin  measured  relative  velocitv  I  and  \  wi  can  estimate  a  and  li  from  the 
dependence  of  measured  currents  on  the  ship  velocity.  We  might  simph  assume  the 
(iirrents  should  he  u  neorrelal  ed  with  the  ship  velocity  and  calculate  a  and  b  from  t  h< 
measured  correlations  and  Eq  !!>  abov<  Even  better,  we  can  ((insider  local  changes  in 
measured  currents  which  ac(ompan\  local  changes  in  ship  velocity.  Define  measures  of 
(hangi  in  current  atui  ship  velocitv  bv 

q  —  I  u  u  I  u  ( 3  ■  1  1  ) 

2 

0  -(V  V  )  V 


where  u  ,u  and  u  art  t  Ik  independent  current  estimates  obtained  steaming  toward, 
steaming  awav  from,  and  while  on  station  respect iveh  .  with  corresponding  definitions  for 
thi  relative  velociiv  estimates  V’  Then  for  each  station  we  can  form  an  estimate  of  the 
((instants  a.h  from  Eq  3>  bv  assuming  that  the  true  values  q  and  Q  are  uncorrelated  ( 
i.e  I  hat  tin  change  in  true  current  is  uncorrelated  with  the  change  in  ship  velocitv  ). 
<ind  bv  approximating  the  true  relative  velocit  v  In  the  measured  relative  velocit  v . 


z'  lq  Q)  IQ  Q) 
(q  Q)  (Q  Q) 


(3.12) 


We  can  then  simph  average  the  estimates  a.  b.  In  practice,  ortlv  estimates  involving 
large  Q  for  which  the  error  signal  should  lie  corresponding!)  large,  were  used,  and  esti¬ 
mates  lying  more  than  .'I  standard  deviations  from  the  mean  were  rejected. 

Table  3. a  shows  the  results  of  this  analysis  for  several  CODE  cruises.  The  cal¬ 
culation  indicates  a  calibration  error  of  about  l'-i.  with  uncorrected  ship  velocity  being 
low.  The  large  value  of  observed  during  CODE  1.  Leg  4  decreases  substantial!)  for 
subsequent  cruises:  this  drop  coincides  with  repair  and  realignment  work  performed  on 
the  gvrorompass  following  Leg  4.  The  data  were  corrected  for  the  calibration  error 
(assumed  to  be  constant  at  1‘7)  arid  for  the  I^eg  4  misalignment  error  before  further 
anal  v  sis 


:i(i 


In  Miimiiiin .  th»  heading  measurement  obtained  from  a  Sperrv -t  v  pt  "trutoin- 
|i<f —  contain^  several  sourees  of  error  of  sufficient  magnitude  to  degrade  estimates  of  tin 
cross-ship  relative  \eioci!\  t  oinjionent .  While  most  of  these  can  be  accounted  for  once 
their  presenti  in  tin  data  is  recognized.  the  rolling  error  arising  from  wave  induced  tilt 
anti  at  1 1  lerat  ions  t  ail  not  - 

3.8.  Calibration  and  Misalignment  Errors 

tin  research  vessel  moves  through  1  lit-  water,  each  component  of  its  horizon¬ 
tal  velocitv  can  varv  over  a  range  from  -  500  cm  set  to  -500  cm  sec.  Since  the  currents 
ttiav  constitute  as  little  a1-  I'd  of  this  signal,  it  is  essential  to  recognize  and  remove  errors 
which  depend  on  the  ship  velocitv. 

>peed  dependent  errors  arise  from  a  wide  variety  of  sources.  Terms  of  order 
|  V  <•)*  in  tin  Doppler  shift .  constant  errors  in  the  oscillator  frequency  or  speed  of 
sounti  used,  splaying  of  the  acoustic  beams,  etc.  have-  the  effect  of  calibration  errors, 
causing  tht'  ship  velocitv  to  be  over  or  underestimated.  Rotation  of  the  transducer  head 
relative  to  the  gyro  reference  line-,  the  misalignment  angle,  causes  the  relative'  velocitv  to 
be  over  or  underrotated  when  referred  to  geographical  axes. 

Suppose  the  data  contain  errors  caused  by  a  calibration  error.  :i.  and  a  misalign¬ 
ment  angle.  S0(,  .  Then  the  measured  horizontal  components  of  relative  velocit  v  will  be 

1  ( 1  -  :i ) l  cosM(  -  ( 1  -  .f) VsirufifV;,  (3.7) 

V  ( I  -  .V)  I  sirnf  0{  -  ( 1  -  rf)  Vcos60( 

where  l  .V  are  the  true  components  of  relative  velocitv.  The  inferred  currents  computed 
from  Eq  1.2  will  thus  be  related  to  the  true  currents  u.v  by 

ti  u  -  bl  -  aV  (3.8) 

V  V  a  l  1>Y 

where 

a  (1  •  Tlsin/-^  (3.9) 

b  (I  -  T)cos <*> /V, ,  I 

If  we  can  use  the  data  to  infer  (a.b).  we  can  determine  the  calibration  and  misalignment 
angle  errors  as 

tan  '(—£—)  (3.10) 

b  •  I 

la*  -  (b-D* 
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rotation  of  the  large  lore  component  of  velocity  into  tin  cros— ship  t omponent .  tin  error 
current'  will  lie  polarized  in  tin  cross-ship  direct  ion.  and  will  increase  with  the  speed  of 
the  ship.  I  nfort  unat  ely .  errors  i'll  of  thi-  order  of  magnitude  cannot  be  ruled  out 

lln  heading  data  for  this  study  were  taken  directly  from  tin  ship’s  Sperr> 
master  gyrocompass.  One  reading  wa-  taken  for  each  acoustic  pulse  I  O.G  set 
sampling  I.  lo  provide  a  heading  measurement,  the  gv  rocompas'  must  sense  its  orienta¬ 
tion  relative  lo  .1  directions.  I  lit  fore  aft  a\i-  of  tin  ship,  the  direction  of  gravitational 
acceleration  and  tin  rotation  a\i'  of  tin  Earth  provide  these  reference  directions.  An 
error  in  the  detection  of  am  one  o|  these  will  yield  an  error  in  the  measured  heading.  A 
complete  discussion  of  the  mechanic'-  of  gyrocompasses  will  not  be  undertaken  here,  and 
interested  readers  art  referred  to  tests  such  as  Arnold  and  Maunder  (19G1)  or  Wrigiev  ft 
»l-  (19(>9).  I  able  T  1  summarizes  some  characteristics  of  measurement  errors  w  hich  are 
present  in  gyrocompasses  of  the  Sperry  type.  The  misalignment  angle  is  the  constant 
angle  between  the  gyrocompass  fore  direction  and  the  fore  direction  defined  by  the 
acoustic  beams,  l  asharn  (197(>)  discusses  a  technique  for  determining  this  angle  which 
recpiires  steaming  a  course  in  a  steady  cross-wind  at  a  variety  of  speeds  and  attributing 
covariance  between  foreward  and  port  ward  relative  velocity  measurements  to  transducer 
misalignment  A  test  run  using  this  procedure  wa-  made  during  CODE:  it  showed  that 
during  the  test  the  acoustic  transducer  was  aligned  to  within  0.2'  of  the  fore  aft  axis  of 
the  ship  as  defined  gyroscopic  ally  by  the  heading.  A  post -processing  method  for  estimat¬ 
ing  both  the  misalignment  angle  and  calibration  error  is  presented  in  the  next  section 
winch  requires  neither  special  test  courses  to  be  run  nor  strong  assumptions  about  the 
steadiness  of  the  wind  and  horizontal  structure  of  the  current  shear  field  to  be  made. 

I  he  latitude  error  and  velocity  error  in  Table  3.1  are  usually  corrected  in  the  course  of 
normal  ship  operations.  The  crew  manually  enters  the  speed  and  latitude  of  the  ship 
into  compensators  located  on  the-  bridge.  Were  this  not  the  case,  these  corrections  could 
have  been  applied  (more  accurately)  in  post-processing.  Acceleration  error,  which  takes 
the  form  of  damped  oscillations  in  the  heading,  at  M  minute  period,  following  a 
maneuver  involving  north  south  accelerations,  can  also  be  corrected  irt  post-processing  if 
the  history  of  the  ship  velocity  is  known  for  a  long  enough  time.  The  wave-induced  rol¬ 
ling  (also  called  quad rantal )  error  depends  on  the  details  of  the  gyrocompass  construc¬ 
tion  arid  I  am  not  aware  of  arty  method  for  removing  it.  The  magnitude  of  0.75r  listed 
for  this  error  comes  from  engineers  at  the  Sperry  Corporation  (Herschel  Porter,  personal 
commuri  icat  ton ) . 
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Table  3.4  :  Heading  Errors 

Error 

Forcing  Form 

Magnitude 

Misalignment  Angle 

installation  error  constant 

') 

Latitude  Error 

A 

y  lanA 

1.6" 

Velocity  Error 

Y  j 

0.8" 

|  RcosA  -  l  | 

Acceleration  Error 

j 

AN  e  *'  sinv 

1  -  0.4" 

[  RQ  "cos^A 

Rolling  Error 

wave  induced  ship  tilt 

- 

0.75" 

where 


R 
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A 

(t.Y) 
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t 


gyrocompass  damping  parameter.  2‘  for  Sperry  Mk  37 
\  g  R  .  the  Schuler  frequency 
radius  of  Earth 
rotation  rate  of  Earth 

IrT 

latitude 

(east. north)  ship  velocity- 

change  in  V  during  maneuver,  assumed  to  occur  quickly  compared  to  It 
time  since  maneuver 


Table  3.4  :  Heading  Errors.  This  table  lists  several  common  errors  which  may  be  present  in  gyro¬ 
compass  data.  Magnitudes  assume  an  operating  latitude  of  38°  N.  a  course  due  north  at  10  knots, 
and  a  maneuver  in  which  the  ship  reverses  course  over  a  time  which  is  short  compared  with  the 
Schuler  period.  NVith  the  exception  of  the  rolling  error,  these  may  be  corrected  either  at  the  gyro¬ 
compass  or  in  post-processing. 
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'I  his  is  idem  ical  to  t  he  res  ults  (-1. M). (-1.22)  since  \  2M  -  1. 

The  two  methods  of  deriving  the  filter  produce  the  same  estimator,  but  each 
provides  a  different  piece  of  supplementary  information.  From  the  optimum  estimator 
calculation  we  obtained  (4.23).  the  expected  mis  error  in  the  filtered  estimate.  On  the 
other  hand,  the  least  squares  formulation  provides  not  only  a  current  estimate  u  (4.34) 
from  the  slope  of  the  drift,  but  also  an  estimate  of  the  intercept  d(..  Since  the  drift  at 
time  zero  must  be  zero,  tl,  is  just  minus  the  error  in  t he  initial  fix  x(1.  By  correcting  the 
initial  position  and  integrating  ( ti  N,)  over  time,  an  improved  estimate  of  the  ship’s 
position  can  be  obtained  This  allows  the  location  to  which  the  estimate  applies  to  be 
determined  more  reliably  than  would  be  possible  by  using  raw  fixe--  alone. 


4.6.  Optimum  Filtering 

The  result  (4.24)  gives  the  expected  rms  error  for  the  estimator  (4.3X).  But 
because  we  have  assumed  the  true  current  is  constant  (or.  more  precisely,  that  4.2a 
holds).  Lq  4. <.  i  predicts  that  this  rrtis  error  will  decrease  to  zero  as  the  filter  length  N  is 
increased.  In  fact,  of  course,  at  sufficiently  large  V  Eq  4.25  no  longer  is  a  valid  assump¬ 
tion.  and  increasing  the  Tiber  length  will  increase  the  error  in  the  estimate  by 
oversmoothing  true  structure  in  the  current  field.  Choice  of  an  appropriate  filter  length 
is  achieved  by  minimizing  the  net  error  due  to  both  sources,  which  can  be  written  as 

(U  u)-  Jsuu( f)(l  W(f))-df  -  f)\V-(f)df  (4.39) 

where  the  first  term  represents  loss  of  signal  due  to  filtering  and  the  second  term 
represents  the  noise  which  passes  through  the  filter.  \V(f)  is  the  transfer  function  of  the 
filter. 

M 

VN  (f)  £  wke  *'*  ika'  (4.4(1) 

I.  M 


which,  after  some  tedious  algebra,  can  be  determined  for  the  optimum  estimator  weights 
(1.22)  to  be 


\N(f) 


(\  ■  2 )sin  \  r  f  At 


Nsin(  N  -  2)tt  f At 


2N(N  -  1 )( N  2) 


sin  r  f At 


(4.4  1 ) 


Since  wo  have  a  model  of  the  noise  spectrum  S  and  can  infer  SUL  from  S-(  and  S(.uflj. 
(■1.39)  can  he  evaluated  as  a  function  of  the  filter  length  N.  and  the  value  of  N  which 
minimizes  ■  (it  u)‘  can  lx-  chosen.  This  technique  also  allows  us  to  test  the  perfor¬ 
mance  of  alternative  filters  by  simply  inserting  tin-  appropriate  transfer  function  into 
(•1.39).  The  results  of  such  a  test  are  shown  in  Fig  4.4  for  3  choices  of  filter  -  the  least 
squares  filter  of  Lq  -1.3-4.  a  low  pass  Tukey  filter  and  a  truncated  si nc  filter.  The  results 
arc  nearly  identical  for  the  three  filters,  with  the  least  squares  estimator  yielding  some¬ 
what  lower  expected  errors  at  short  filter  lengths  and  achieving  minimum  error  more 
rapidly  than  the  other  low  pass  filters,  and  the  other  filters  removing  somewhat  less  sig¬ 
nal  in  the  case  of  overfill ering.  In  all  cases,  the  mean  square  error  drops  steeply  to  the 
minimum,  then  rises  slowly  as  the  filter  length  is  increased  past  t fie  optimum  value. 
Since  we  don't  know  with  great  confidence  where  the  signal  truly  begins,  the  data  were 
processed  conservatively  by  using  a  filter  of  length  V- 30.  Fig.  -4 .-4  suggests  an  rms 
uncertainty  of  2-3  cm  sec  in  \.  t  fie-  alongshore  current,  and  4-5  cm  sec  in  u.  the  cross- 
shore  current,  should  be  expected. 

4.7.  Vertical  Current  Profiles 

The  analysis  to  this  point  has  focussed  on  predicting  the  current  u(zk)  at  a  sin¬ 
gle  depth  zk  from  acoustic  measurements  for  that  depth  and  from  navigational  data.  To 
estimate-  the  current  at  another  depth,  of  course,  the  filtering  procedure  could  simply  be 
repeated  using  the  acoustic  data  from  the  new  depth.  However,  w'e  have  seen  that  the 
acoustic  data  is  much  less  noisy  than  the  navigation  data,  so  that  while  30  minutes  of 
combined  measurements  are-  required  to  determine  the  absolute  current  at  any  level, 
vertical  shears  between  levels  can  be  determined  from  5  minutes  of  acoustic  data  alone. 
Estimating  the  current  at  one  depth  using  30  minutes  of  data,  the  vertical  profile  of 
currents  could  then  be  determined  from  the  shear  over  the  shorter  5  minute  period, 
retaining  shorter  scale  variability  in  the  current  profiles.  To  obtain  the  current  profiles 
for  this  study.  30  minute  time  series  of  vertically  averaged  V  were  used  to  estimate  the 
average  current  over  a  vertical  slab  of  ocean.  The  vertical  deviations  from  this  average 
were  then  determined  from  the  central  5  minute  average  profile  of  V.  I  sc  of  this  slab 
determination  results  in  ecpial  smoothing  of  currents  at  all  levels  by  the  30  minute  filter. 

4.8.  Summary 

Flic  goal  of  this  section  was  to  use  the  imperfect  measurements  of  position  and 
relative  velocity  described  in  Chapter  3  to  best  extract  the  current.  Our  approach  has 
been  to  use  each  pair  of  adjacent  fixes  and  the  average  relative  velocity  between  them  to 
produce  a  single  noisy  measurement  of  the  current.  We  have  shown  that 

(a)  the  spectrum  of  such  measurements  shows  a  strong  spectral  gap  separating  two 


peak*. 

Ui)  tin  high  frequency  peak  is  rimclt-l let)  well  over  most  of  its  range  by  assuming  it 
arises  from  white  noise  error  in  the  position  measurements. 

|c]  the  nois<  in  an  estimate  it  formed  by  averaging  N  adjacent  measurements 
decreases  as  I  N.  anti 

(dl  ati  estimator  formed  by  increasing  t  tie  time  between  fixes  is  equivalent  to  (c). 

By  assuming  the  true  current  obeyed  ( 4.2”*)  over  the  filtering  interval,  we 

(e)  derived  t  lie  estimator  (-1.11  or  -1.22)  w  hich  gives  the  minimum  mean  square  error 
(1.2-1 ).  and 

(f)  showed  that  this  estimator  is  equivalent  to  making  a  least  squares  estimate  of  the 
drift  of  a  water  parcel. 

By  assuming  the  true  current  had  variability  described  by  the  measured  spectrum  after 
removing  the  spectrum  of  white  noise  fix  errors,  we 

(g)  tested  the  robustness  of  the  optimum  estimator,  and  found  it  to  be  not  much 
better  or  worse  titan  other  possible  low  pass  filters,  and 

(it)  found  the  dependence  of  ruts  error  in  the  filtered  signal  on  filter  length. 

We  then  commented  on  how  to  extend  the  single  depth  analysis  to  vertical  profiles. 


COMPARISON  W  I  TH  MOOR  PI)  CCRRLNT  ML'TKRS 


5.1.  Introduction 

The  COI) I.  2  moored  array  of  current  meters  provides  an  independent  set  of 
measurements  taken  while  the  Doppler  acoustic  log  was  being  operated.  By  comparing 
the  current  meter  measurements  with  those  inferred  from  tin-  shipboard  Doppler  as  it 
moved  through  the  current  meter  array  we  may  learn  about  the  relia  l>i  lit  \  of  the  ship¬ 
board  Doppler  technique.  Furthermore,  the  fixed  current  meter  movable  DAL  geometry 
y  ields  simultaneous  measurements  along  a  continuum  of  spatial  separations.  This  means 
that  s  pat  ia  I  scales  cart  be  investigated  directly  with  the  combined  data  set.  without  tin 
inherent  mixing  of  space  and  time  variability  to  which  survey  measurements  alone  art 
subject . 

The  CODL-2  moored  array  (Fig.  5.1)  was  deployed  in  March  1982  by  scientists 
from  Scripps  institution  of  Oceanography  and  Woods  Hole  Oceanographic  Institution, 
and  was  recovered  in  July-August  1982.  Moorings  were  set  along  4  lines  (designated  1. 
V  C.  and  R  from  north  to  south)  perpendicular  to  the  coast,  in  water  depths  of  60m 
(except  the  I  line).  90m.  l.'lOm  and  400rti  (C  line  only).  Details  of  mooring  locations  and 
instrument  depths  are  given  in  Table  5.1. 

Data  from  these  moorings  were  kindly  made  available  by  R.  Davis.  C.  W  inant 
and  R.  Beardsley  as  vector  averaged  hourly  currents. 

5.2.  Comparison  within  1  km 

Direct  comparison  between  the  DAL  and  current  meter  measurements  was  made 
a^  follows: 

Whenever  the  ship  location  (  from  LORAN-C  positions  corrected  as  described  in 
Chapters  5  and  4  )  was  within  1  km  of  the  nominal  location  of  a  mooring  (  also 
corrected  for  constant  LORAN-C  offset  ).  the  DAL  profile  was  linearly  interpolated  to 
tin  current  meter  depths,  and  the  hourly  current  meter  values  were  linearly  interpolated 
to  tin  DAL  measurement  time.  Current  meters  which  were  either  too  deep  or  too  shal¬ 
low  for  valid  DAI.  estimates  were  not  used.  Comparisons  are  made  using  the  alongshore 
and  cross-shore  components,  where  the  alongshore  direction  is  taken  to  be  .'t  1 71  every¬ 
where.  except  for  the  I  moorings,  where  it  is  taken  to  be  ()"T. 


Fig.  5.1  :  CODE  2  (1982)  moor 
instrumentation  at  each  mooring 
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current  meter  array.  Table  5.1  provides  details  on 
Also  shown  is  the  location  of  wind  buoy  NDBO 
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T able  5 

1  CODE  1 

MOORED  CURRENT  ARRAY 
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Table  5.1  :  CODE  2  moored  array.  S  (SS)  denotes  a  current  meter  string  which  is 
surface  (sub-surface)  moored.  \  ACM  (VMCM)  denotes  a  Vector  Averaging  (Measuring) 
Current  Meter. 


Comparisons  were  performed  for  the  entire  CODL  2  data  set  over  a  wide  range 
of  processing  options.  Filter  lengths  of  half  hour,  hour  and  two  hours,  filters  of  th<  least 
squares.  Tukey  and  endpoint  type,  siah  and  single  level  determinations,  as  well  as  \ari- 
ous  calibration  corrections  were  all  tested.  No  alternatives  to  the  half  hour,  least  squares 
filtered,  slab  averaged,  calibration  corrected  processing  consistently  yielded  improvement 
in  tin1  comparison  statistics.  A  substantial  improvement  in  comparison  results  was 
obtained,  however.  In  an  ad  hoc  change  from  the  nominal  depths  of  the  1)AL  range  bins. 
By  treating  the  acoustic  data  in  range  bin  n  ]  as  if  they  had  come  from  range  bin  n.  the 
rrtts  difference  between  current  meter  and  I)Al.  measured  currents  was  dramatically 
redu<  •  d  (typically  3(l(<  reduction  in  the-  variance-  of  the-  difference-  signal).  Subsequent 
measurements  by  Regier  (personal  communication)  confirmed  the-  presence  of  hardware 
sources  fe>r  such  delays  in  the-  Arnetek  ele-ctronics.  The  re-sults  presented  here  use  the 
corrected  depths  for  the  measurements. 

'I  he-  i it t  ercomparison  results  are-  shown  graphically  by  the  scattei  plots  in 
I  ig  ‘>.2:  a  statistical  summary  of  the  results  is  given  in  Table  5.2.  In  the  table  there  are 
two  columns  which  count  the-  number  of  comparison  points:  the  first.  Nlot.  gives  t he¬ 
lot  a  I  number  of  comparison  pairs  (one-  DAL  me-asurerrient  per  5  minutes)  and 
corresponds  to  the-  number  of  points  shown  in  the  scatter  plots.  while  the  secono. 
labelled  NIIld  .  cetunt s  only  those  samples  for  which  distinct  pairs  of  hourly  current 
meter  averages  were-  used  in  deriving  a  time-interpolated  current  to  compare-  with  the 
DAI.,  and  is  the-  metre  applicable  figure  for  determining  significance  levels. 

The  rtte-an  cross-shore-  current  component  u  at  each  depth  as  measured  by  the 
DAL  differs  by  ().(i  cm  sec  or  less  from  the  mean  over  simultaneous  measure-merits  by  the 
metore-d  current  meters:  however  the  means  for  the  alongshore  component  \  differ  by  up 
to  1  .(>  cm  sec.  with  the  DAL  consistently  showing  stronger  dowricoast  flow  than  the 
enrre-nt  meters.  There-  is  a  strong  cross-shore  and  vertical  gradient  in  the  mean 
alongshore-  component  of  current  over  the  shelf,  as  we-  shall  see  later,  and  the  differences 
in  \  may  arise  freiru  small  e-rreirs  in  our  knowledge  of  the  positions  for  the  measurements 
be-ing  e  ompared.  (  -e»r re-lat  iems  be-twe-en  the  two  sets  of  measure-merits  are  high,  ranging 
from  *2  for  u  anel  .95-. 97  for  v.  Standard  deviations  of  the-  differences  are  -l.l-5.-t 
etn  see  in  ii  and  1  t  cm  see  in  \:  recall  that  errors  of  4-5  cm  sec  in  u  and  2-3  cm  see 
Hi  \  were  expected  in  the  DAL  currents  from  the  considerations  in  Chapter  -1.  The  size 
of  the  differences  between  DM.  and  moored  current  meter  measurements  of  currents  was 
found  to  increase  bv  small  but  statistically  significant  amounts  as  either  the  ship's  speed 
or  distanee  from  the  mooring  increased,  indicating  that  both  measure-men t  errors  and 
geophysical  variability  contribute  to  observed  differences.  No  tendency  was  found  for 
the  mean  or  fluctuating  differences  to  have  a  preferred  orientation  relative  to  the  ship. 
The  variance-  of  the  DAL  measurements  is  consistently  larger  than  that  of  the  current 
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Tame  "i.'J  :  Current  tneter  UAL  comparison  statistics.  Whenever  the  DAL  was  within  1  km  of 
,1  isi'iirirt'p  comparison  of  DAL  and  moored  current  meter  measurements  was  made.  Averages  and 
.  ur’.ahre^  to-  each  instrument  are  shown,  as  are  covariance,  correlation  and  standard  dev  iation  of 
•he  difference  between  measurements  from  each  instrument.  Also  shown  is  a  comparison  between 
•  ne  or.iv  ilosrlv  spaced  pair  of  moored  instruments  in  (  UDK.  a  \  ACM  and  a  \  MC  M  both  at  10tn 
;n  .ir  t  '■ 
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Fig.  5.2a  :  Comparison  of  all  DAL  current  measurements  (alongshore  component) 
made  within  1  km  of  a  current  meter  mooring  with  the  simultaneous  current  meter 
measurement.  At  both  20m  and  35m  depth,  over  1000  pairs  of  measurements  are 
plotted.  See  Table  5.2  for  a  statistical  summary  of  these  results. 
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Iii  asse-ssing  t  h  i'-  comparison  it  is  useful  lo  know  how  two  "standards"  compare 
with  one  another.  Therefore  direct  comparison  was  also  made  between  the  only  pair  of 
current  meters  in  ('OI)I  2  which  were-  located  at  the  same- depth  and  closeh  spaced  hor- 
i/ontally.  one  a  \  A(\l.  the  other  a  \  \1('M.  Both  were  located  1  On i  below  surface 
moorings  in  approximately  9(lm  of  water  along  the-  ('  line:  the-  separation  between  moor¬ 
ing'-  wa-  approximat cl\  lOftm.  considerably  les*-  than  the  typical  distance'  the  ship  occu¬ 
pied  during  comparison  times.  The-  comparison  was  made  for  the  common  time  period  of 
121.5  days  beginning  at  1500  (I  Ml  on  21  March  I9S2.  ’The*  result^  are  also  presented  in 
I  able  5.2.  The  differences  between  the  measurements  are  surprising!)  large  -  mean 
differences  of  ().“»  and  2. a  cm  see  and  standard  deviations  of  the  differences  of  4.2  and  5 . 3 
cm  see  in  u  and  \  respect  i\ el\ .  These  differences  do  not  correlate  well  with  the  currents, 
and  sei  are  not  cas'd)  explained  in  terms  of  fouling,  broken  fans,  calibration  errors,  etc. 

The  comparison  results  are'  quite  encouraging  and  demonstrate  that  meaningful 
current  measurements  cari  be  obtained  from  the-  shipboard  I)AL.  More  aecurate  rneas- 
ure-me'tits  of  position  and  heading  should  reduce  errors  further.  An  important  obstacle' 
which  is  more-  difficult  to  surmount  is  asse-ssing  and  eorre’Ct ing  for  variability  in  the- 
acoustic  environment  below  the  ship.  especially  the1  influence  that  bubbles  have  on  sound 
speed  and  thus  on  the'  re-lat  ion  between  l>opp  ler  shift  and  relative  velocit)  (cf  Eq.  2.1). 

5.3.  Comparison  at  Larger  Separations 

Moored  current  meters  continuously  measure'  the-  current  at  a  single  location. 
The  shipboard  DAI.  measures  current  at  the  ship's  location.  Thus  as  the  ship  moves  in 
the-  vicinity  of  a  curre-nt  rnele-r  the  two  instrument'-  form  an  adjustable  array  and  allow 
simultaneously  measured  current'-  to  be-  compared  at  a  variety  of  spatial  lags.  During 
CODE  the  C '-line  was  the  region  fe»r  which  shipboard  sampling  was  most  intensive 
(Appendix),  and  thus  statistical  reliability  highest.  Eig.  5.3a  shows  how  the  mean 
currents  measured  by  the  DAL  varied  across  the*  shedf  along  the1  C  line.  The  origin  is 
taken  at  the  <'3  mooring.  A  strongly  varying  mean  field.  as  observer!  for  v(x).  increase’s 
the  apparent  noise  for  intercomparison  experiments  such  as  the  one  just  discussed. 
Divergence  in  the  me-an  cretss-slieire-  current  u.  e  xpected  in  the  presence  of  ceiastal  upwed- 
ling.  is  clearly  seen.  Fig.  5.3b  shows  how  the-  correlation  between  DAL  and  moored 
measurements  at  03  fall  off  as  a  function  of  offshore-  separation  of  the-  ship  from  the 
mooring.  For  u.  correlation  between  the  instruments  falls  to  values  not  significantly  dif¬ 
ferent  from  zero  at  separations  less  than  the'  mooring  separations,  while  for  v.  the  cross- 
shore-  ceirre'lat  ion  length  is  substantially  longer.  Fig  5.3e  shows  that  the'  rttis  differe-nce 
bet w e-en  currents  grows  rapidly  as  a  function  of  cross-shelf  separation. 
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fig.  5  3b  .  Correlat ion  between  simultaneous  DAL  and  C3  current  meter  currents  as 
a  function  of  cross-shelf  seperation  from  the  C3  mooring.  Dashed  line  gives  95(’t 
confidence  level  for  correlation  over  available  samples.  Stars  show  C3  correlation  with 
current  meters  at  C2  and  C-1. 


Fig,  6.2b  :  Secondary  grid 
100  fathom  isobaths,  and  ; 
Appendix  for  times  occupied. 


Fig.  6.2a  :  Primary  grid  of  CTD  stations  for  CODE.  Stations  on  each  line  are 
numbered  consecutively  from  station  nearest  shore.  The  Irish  Gulch  line  was  sampled 
by  CTD  only  during  1982.  The  1982  Central  line  was  some  3  km  south  of  the  1981  line. 
See  Appendix  for  times  occupied. 
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Fig  6.1  :  Alongshore  component  of  surface  wind  stress  (dynes/cm2)  at  NDBO  46013. 
computed  from  hourly  wind  measurements.  The  hourly  stress  values  were  low  passed 
using  a  filter  with  a  40  hour  half  power  point. 
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I  ig.  6.6  i-  a  map  show  in"  <1  principal  axis  analysis  of  current  fluct  nations  about 
tin  mean  at  each  ('ll)  station  lor  a  depth  of  30m.  A-  noted  by  Kuridu  and  Allen 
(I97(il  I  lit  re  is  a  si  rong  alignment  of  the  principal  axes  near  the  coast,  while  farther  from 
short  tin  fluctuation-  approach  isotropv  and  so  t  lie  principal  axes  lit'  along  random 
direction-  (her  tin  shell,  t  ht  principal  axis  direction  rotates  offshore  in  the  same  sense 
a-,  but  somewhat  more  strongly  than,  t  lit-  isobaths.  Along  the  Arena  line,  representing 
tin  flow  around  a  turner,  this  rotation  is  especiallv  pronounced  and  fluctuations  quickly 
become  si  congest  in  a  direction  which  is  distinct  ly  away  from  the  coast. 

6.3.  Discussion  and  Summary 

"Iht  shipboard  current  measurements  provide  a  view  of  the  mean  circulation 
w  hicli  is  very  well  sampled  in  space,  but  not  continuously  sampled  in  t  ime.  It  is 
apparent  from  Fig.  6.)  that,  especially  during  19N2.  sampling  may  be  biased  toward 
periods  of  low  winds.  The  bia-  in  the  mean  currents  due  to  shipboard  sampling  of  the 
wind  forcing  was  estimated  at  each  CTI)  station  from  a  simple  regression  analysis  of  the 
moored  current  measurements  on  the  wind  stress  at  NDBO  46013.  The  predicted  bias 
was  small  and  did  not  significantly  affect  the  structure  of  the  fields  presented  above. 

The  dens<  spatial  sampling  provided  by  the  shipboard  DAL  clearly  delineates 
tin  extent  and  intensity  of  features  in  the  coastal  upwelling  circulation.  As  expected, 
tin  mean  cross-shelf  circulation  shows  an  Lkman  layer  of  offshore  transport  near  the  sur¬ 
face  with  return  flow  below.  Out  to  the  shelf  break,  the  layer  of  offshore  flow  deepens 
with  distance  from  the  coast.  The  poleward  undercurrent  is  clearly  seen  in  the  mean 
alongshore  current:  it  appears  strongest  near  the  shelf  break,  and  surfaces  near  the  roast. 
A  near  surface  equatorward  jet  is  apparent  in  the  mean  alongshore  current  across  each 
hydrographic  line,  but  since  it-  core  moves  offshore  from  north  to  south,  it  is  artificially 
smoothed  out  when  means  are  calculated  a-  a  function  of  offshore  distance  alone.  The 
current-  show  significant  divergence  in  both  alongshore  and  cross-shore  components, 
from  which  the  mean  vertical  current  can  be  calculated.  Mean  upwelling  appears  to 
occur  over  the  entire  shelf,  providing  enough  vertical  transport  to  balance  the  predicted 
mean  Lkman  transport  in  the-  (’OI)L  region.  Fluctuations  about  the  mean  current  are 
strongly  polarized  in  the  alongshore  direction  near  the-  coast  and  essentially  isotropic  far 
Irom  tin1  ciia-l  However,  flow  near  the  coastal  corner  at  Ft.  Arena  shows  fluctuations 
polarized  away  from  the  coast.  Ferhap-  t  he  most  interesting  finding  of  the  CODL  DAL 
measurements,  however,  is  that  the  relative  simplicitv  of  the  average  field  does  not  carry 
over  to  the  svnoptu  field,  a-  the  mapping  surveys  discussed  in  the  next  chapter  show. 


lor  t  Ik  dv  nami<  height  (e.g.  Reid  and  Manly  la.  1970)  which  an  required  lo  infer  tin 
gcostrophic  current  in  shallow  water  An  example  of  the  limitation1-  of  such  scheme--  i- 
given  in  sect  ion  7.1. 

Significant  divergence  du  dx  is  observed  in  the  mean  cross-shore  current  over 
the  shelf  (Fig.  b.ilf).  positive  in  tin  upper  water  column  and  negative  deeper  down.  In 
tin  upper  water  column  over  the  shelf,  the  cross-shore  divergence  averages  .'5  ■  1(1  1  set 
thi-  agrees  with  Davis'  (19Mb!  estimate  of  surface-  divergence  from  drifter  measure¬ 
ments.  If  t  lit  divergence  in  tin  alongshore-  current  can  be  neglected,  mass  conservation 
and  tin  data  in  Fig.  (>..‘{f  imply  that  tin  mean  vertical  (upwellingl  current  will  be 
directed  toward  the  surface,  with  a  maximum  along  the  du  dx  -  (I  contour.  Although 
tin  I)-\l.  data  oo  not  reach  the  surface  or  the  bottom,  integration  over  the  data  shown 
suggest-  a  mean  upwelling  rate  at  the  midwater  maximum  of  at  least  0.7-  10  *  cm  sec. 
or  (>  m  dav.  averaged  over  the  inner  15  km  of  shelf.  The  implied  vertical  mass  transport 
over  this  part  of  the  shelf  would  be  about  10'  m*  dav.  an  amount  equal  to  the  offshore 
likinai:  transport  forced  In  a  mean  surface  wind  stress  of  1  dyne  crir. 

Significant  divergence  in  u  is  also  indicated  in  the  vicinity  of  tin-  deep  maximum 
in  u  near  the  shelf  break:  whether  this  result  represents  true  physics  or  simply  an  over- 
optimistic  estimate  of  statistical  reliability  is  not  known.  The  limited  evidence  from  the 
moored  instruments  at  C5  do  show  an  increase  in  u  in  this  region. 

Averaging  over  y  in  the  above  analysis  implicit )y  assumes  that  alongshore  mho- 
mogeneitv  in  the  statistics  is  unimportant.  This  i-  done  in  order  to  provide  increased 
densitv  of  observations  in  a  broad  overall  section.  Aspects  of  alongshore  variability  in 
the  statistics  can  be  examined  by  comparing  the  average  sections  obtained  along  each 
(  II)  line.  The  sections  are  shown  in  figs,  (i.-la.c  arid  the  uncertainty,  significant ly 
increased  over  that  of  the  y  -averaged  section,  is  shown  in  Fig.  (>.4b.  I  he  structure  of  the 
mean  alongshore  current  v  i-  considerably  more  jet-like  in  these  individual  sections  than 
in  the  overall  average  Fig.  ( j . '{ b .  f  rom  north  to  south,  the  tore  of  the  jet  accelerates  and 
migrates  offshore.  The  latter  result  is  consistent  with  the  fact  that  the  strongest  offshore 
flow  appears  in  Fig.  b.-la  to  he  associated  with  the  core  of  the  jet. 

The  alongshore  structure  in  the  mean  Held  implies  a  non-zero  dv  dy.  negative 
onshore  (positive  offshore)  of  the  core  of  the  jet.  which  contributes  to  the  total  diver- 
genie  responsible  for  the  mean  upwelling.  Because  the  core  of  the  jet  migrates  offshore, 
dv  dv  is  a  function  of  the  alongshore  coordinate  y.  Rather  than  attempt  a  simple  y 
independent  representation  of  dv  dy  to  combine  with  du  dx  from  Fig.  G-.'lf,  dv  dy  w  a- 
estimated  from  two  pairs  of  CTI)  lines  and  the  results  shown  in  Fig.  (>.5.  The  effect  of 
dv  dv  is  to  substantially  increase  the  inferred  maximum  w  over  the  shelf  and  to  decrease 
1 1  bey  on d  t  he  shelf  break. 


Hi.  mean  alongshore  and  cross-shore  currents  art  shown  in  I  ig.  dd'.li.  \lanv  of 
tin  classical  feature-  expected  for  a  coastal  npwelling  flow  field  (  e.g.  Allen  I  19*0). 
Iluver  (  1  )  )  can  be  seen  dearlv  in  the  average  fields.  (  (insistent  with  the  observed 
cquatorward  mean  wind  stress,  a  near  surface  laser  of  offshore  transport  is  evident  at  all 
off-horc  locations  in  fig.  (i.lib.  The  thickness  of  tiiis  laver  increases  from  les-  than  20m 
near  tin  coast  to  about  50m  near  the  shelf  break.  Weak  onshore  flow  occurs  every  when 
below  this  later,  except  in  a  region  over  the  shelf  near  t  tie  bottom  of  the  I)A1.  resolution 
I  where  the  difference  from  zero  is  generallv  not  significant  according  to  fig.  (did  ).  It  i- 
thi-  pattern  of  t  ro---shore  circulation  which  produces  the  npwelling  of  cole  water  near 
the  (mast.  Near  tin  surface,  tin-  alongshore  mean  flow  is  characterized  b>  an  equator- 
ward  jet  extending  from  midshelf  to  beyond  the  shelf  break,  flowing  counter  to  the 
cquatorward  mean  wind  stress,  a  poleward  undercurrent,  strongest  near  the  -tielf  break, 
surfaces  near  t  he  coast . 

I  luct  uat  ions  in  it  and  v  about  their  mean  values  (  fig.  G.lic  )  an  at  least  a- 
largt  a-  the  mean-  themselves,  far  from  shore,  tin  fluctuations  in  u  and  v  are  equal!) 
energetic.  As  tin  coast  is  approached,  the  fluctuations  in  u  arc-  damped  dramatical!), 
fliese  results  echo  the  schematic  picture  of  the  npwelling  flow  field  inferred  oxer  the  past 
two  decade-  from  point  current  measurements  and  hydrograph)  (fig.  1.1  I. 

Cross-shore  and  vertical  gradients  in  the  mean  field-  were  computed  it)  smooth¬ 
ing  first  difference-  between  adjacent  (x.z)  bins  with  a  (.'>•.'{)  triangular  weight  filter. 
The  results  are  shown  in  fig-  die  and  f.  In  the  interior.  awav  from  surface  and  bottom 
stresses,  t  tie  vertical  shear  in  n  is  small  (generall)  less  than  10  “  sec  1  |.  Above  40m.  u 
becomes  increasing!)  sheared  as  the  surface  is  approached,  as  expected  in  the  presence  of 
wind  forcing,  for  v.  bv  contrast,  substantial  vertical  shear  is  evident  over  a  much 
greater  dept  it  range  ,  as  expec  ted  if  the  upvvclled  densit)  surfaces  are  geost rophical!)  bal¬ 
anced.  A  rough  est  imate  of  t  tie  mean  geost rophic  shear  can  be  made  from  Miner's 
report  of  the  average  of  17  post -t  rans'u  ion  In  drographic  sections  along  the  Central  line 
during  I9M  (Iluver  I9MI  Substituting  measured  isopycnal  slope's  and  iSp  dz  from 
Miner's  fig.  1  into  t  he  thermal  wind  equation 

civ  g  c)  p  dz 

ciz  f,  f  dz  dx 

)  ield-  a  mean  siiear  of  1.5  ■  10  "  a-  a  broad  average  for  deptiis  greater  than  50m  over 

the  shell,  'flu-  is  in  rough  agreement  with,  though  slightlv  smaller  than,  the  direct  1) 
measured  values  A  detailed  comparison  over  paired  DAL  CTI)  measurements  is  indi¬ 
cated.  and  will  be  undertaken,  to  find  the  small  deviations  from  geost  rophy  which  are  of 
most  dynamical  interest.  Note  t  fiat  such  direct  comparison  of  measured  and  geostTophir 
shears  obviates  the-  need  to  invoke  an)  of  the  neccssaril)  (id  hoc  extrapolation  schemes 
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6.1.  Data 

During  1  in  -'priii"  and  summer  upwelling  seasons  of  198]  and  1982.  in  conjunc¬ 
tion  with  ('Ol)l..  some  71  da\s  of  enrreni  data  were  collected  in  the  coastal  and  offshore 
water*  of  Northern  California  using  the  Doppler  Acoustic  Log  on  hoard  the  K  \ 
Weronia.  The  data  were  collected,  for  the  most  part,  during  normal  ships  operation 
relating  to  other  aspects  of  CODE.  primarily  during  the  hydrographic  surveying  cruises. 
Figure  (>.  |  shows  the  periods  of  Doppler  log  operation  in  the  CODE  area,  together  with 
the  alongshore  component  of  the  surface  wind  stress,  computed  from  hourly  wind  meas¬ 
urement'-  at  Coast  Cuard  Buoy  .NT) BO  TGOl.'i  and  -10  hr  low-pass  filtered  (this  data  was 
kindly  made  available-  by  Oeorge  Halliwell).  Figures  (i. 2  show  the  CTD  stations  which 
comprised  the  primary  and  secondary  survey  grids  during  Doppler  Log  i  ’ration.  I  he 
sampling  history  at  each  CTD  station  is  shown  in  the  Appendix.  Some  dditional  data 
wen  obtained  in  1982  during  mooring  deployment  and  recovery  cruises.  In  order  to 
focus  on  coastal  dynamics  during  the  upwelling  season,  data  obtained  prior  to  the  1982 
spring  transition  (Legs  2-’>)  were  not  included  in  the  analysis. 

We  shall  first  discuss  the  averaged  measurements,  then  consider  spatial  variabil¬ 
ity  in  the  averages,  arid  ((include  in  the  next  section  with  ari  examination  of  the  synoptic 
fields  from  w  hich  t  lie  average*  are  formed.  The  notat  ion  w  ill  be  standard,  with  x.y.z.t 
indicating  onshore,  upcoast.  vertical  (positive  upwards)  and  time  coordinates  respec- 
tivelv.  Currents  are  resolved  into  cross-shore  and  alongshore  components  (u.v)  along 
each  primary  CTD  line,  with  the-  alongshore  direction  being  defined  as  l)  1  for  the  Elk 
and  Irish  Culch  lines.  'I  for  the  Arena  line,  and  ,''17f  T  for  the  North.  Central  and 
Boss  lines. 

6.2.  Mean  Fields 

Figures  6.M  present  statistics  of  the  coastal  current  measurements  as  functions  of 
offshore  distance  and  depth.  The  results  were  obtained  by  first  averaging  the  data  in 
spare-time  bins  of  size  (dx.dy.dz.dt)  (2. a  km.  a  km.  (>.,>  m.  1  day)  centered  along  each 
primarv  CTD  line  in  Fig.  6.2a.  then  averaging  over  y  and  t.  The  number  of  observation 
bin*  at  each  (x.z)  is  shown  in  Fig.  6.3a.  Each  bin  has  been  treated  as  an  independent 
sample  (  the  short  lime  scale  was  chosen  because,  in  general,  resampling  of  a  station 
after  more  than  I  but  less  than  2  or  days  was  not  done  at  random,  but  by  design  when 
w  ind  conditions  had  shifted  significant  ly  ). 
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of  the  mean  for  Fig.  6.4a. 


6  4b  Standard  error 
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7.1.  Introduction 

Tin  prec  ed  iii"  desrript  ion  of  tin  average  current  field  is  deceptive  in  it-  simpli- 
rit  \  The  average  field  i-  smooth.  lluct  uat  ions  about  the  local  means,  however,  art'  at 
least  a1-  large  as  the  mean-  themselves  (fig-.  (>..”>).  "1  he  average  field  i-  then  not  neces¬ 

sarily  th('  tvpical  field.  In  1  hi-  chapter  we  shall  present  maps  of  currents  measured  dur¬ 
ing  individual  DAI.  surveys.  These  maps  reveal  a  high  degree  of  spatial  variability  in 
the  t  \  pica  I  field.  Within  this  gencrallv  complex  flow  field,  some  features  appear  to  recur. 
Two  of  these,  the  relaxation  from  strong  wind  forcing  and  the  tongues  of  cold  water 
injected  from  coastal  into  offshore  waters,  will  be  discussed  briefly. 


7.2.  Synoptic  it  y 

To  what  extent  will  maps  made  from  the  (01)1.  Doppler  data  be  synoptic.' 
Davis  (19Mb)  has  estimated  a  Lagrangian  time  scale  of  1.5  days  for  current  fluctuations 
from  surface  drifter  measurements,  and  points  out  that  this  is  significantly  shorter  than 
tin  a  day  Lulerian  time  scale  determined  from  24-hour  filtered  current  meter  measure¬ 
ments.  indicating  the  presence  o  strong,  relatively  long-lived  structure  in  the  current 
field.  It  should  be  noted  that,  since  the  ship  travels  through  spatial  structure  an  order 
of  magnitude  faster  than  do  the  drifters,  the  I)AL  time  scale  will  be  much  shorter  than 
e\  en  t  he  drifter  t  true  scale. 

To  determine  the  time  scales  of  the  energetic  patterns  of  current  variability,  spa¬ 
tial  empirical  orthogonal  functions  (KOFs)  were  determined  from  the  CODL-2  hourly 
moored  current  meter  data  from  2am  depth  and  the  decorrelation  times  (Davis.  197(>) 
calculated  for  the  time  varying  amplitudes  of  each  LOF.  As  expected  the  more  energetic 
modes  have  larger  spatial  scales,  while  the  les-  energetic  modes  contain  the*  shorter  scale 
structures.  Decorrelation  times  for  the  energetic  modes  were  also  longer.  The  first  four 
KOI  s  explain  79'  i  of  the  total  variance  (4*,'<.  I9'i.  7'j.  and  (i'c  respectively)  and  have 
an  energy -w  eight  ed  decorrelat  ion  t  ime  scale  of  4>  days  (3.5.  9.5.  2. a  and  2.4  day  s  respec  ¬ 
tively).  I  or  the  coastal  DAL  surveys  to  be  shown,  the  average  time  separating  pairs  oi 
measurements.  t  t  ..  .  ranges  from  less  than  1  f*  to  about  1  4  of  this  decorrela¬ 

tion  time.  Thus  although  the  DAL  surveys  are  not  a  snapshot,  the  features  seen  arc 
much  more  indicative-  of  spatial  structure  than  temporal  variability  . 


7.3.  DAL  Atlas 


I  r .  1  ig-  7.1  to  7.1  I  we  present  an  at  la-  of  t  hr  phenomena  observed  in  t  lie  coastal 
current  field  during  tin  - hi |>lioar<l  mapping  effort  of  CODL.  The  current  measurements 
c  rnc-  Iron,  a  cieptii  ol  2*mh.  for  clarity.  measurement'-  separated  by  less  than  ;i 
minimum  di-tami  Itvpieallv  3  kmt  haw  beet,  averaged.  I  lie  currents  have  been  over¬ 
laid  on  salelbti  IH  itnajc'  of  s<-,;  -airfare  temperature  whenever  relatively  cloud-free 
image-.  <  los<  enough  in  time  to  the  ship  surveys.  were  available  from  the  >cripp'  Remote 
'sensing  I  aeilitv.  I.aeh  image  cotitain-  a  le-gcnel  showing  the  date-  of  both  the-  satellite 
linage  and  s|,ip  survev.  Tie  mark'  have  been  placed  evert  half  degree  of  latitude  aiiei 
longitude  anil  .i  -e  aie  vector  for  the  currents  i-  shown.  The-  length  of  ant  e  urrent  arrow 
e  orrespoinb  to  tin  elisplae  e-men  t  lit  that  current  over  a  1  .(>  hour  period  (9.3  hours  tor  t  he 
2  large-  survev -j  I  i g  (i.'J  eaii  be  overlaul  on  each  image  to  relate  the  feature's  seen  to  the 
li.it  tiv  me  t  rv  I  tg  (ol  can  be  used  to  gauge-  the  wind  forcing  during  each  survev.  The 
hydrographic  measurements  taken  during  each  survev  are  described  in  a  series  of  Data 
Re  port'  bv  liny  er  ami  co-worker-  IOSI  Reference  series). 

I  he  satellite  IR  shows  that  sea  -urlaet  temperature  varies  strong!)  in  both  the- 
aiong-hore  and  e  Toss-shore  direct  ion.  Although  such  alongshore  variability  has  long  been 
tinted  iti  satellite  images  (  Derust ein .  e/  <il.  19721.  its  relevance  as  an  indicator  of  coastal 
i  ire  illation  and  exchange  processes  has  remained  in  doubt  Without  in  situ  measure¬ 
ment'  tin  depth  extent  of  the  feature'  and  strength  of  the  associated  circulation  was 
line  It  at 

flu  l)\i  current  measurements  exhibit  a  strong  v  isual  correspondence  with  the 
satellite  derivee!  surface  temperature  field.  Where  this  is  true,  the-  features  in  the-  satel¬ 
lite  image-'  lannot  -imply  be-  dismissed  as  skit:  effects  on  the  sea  surface.  Rather,  the 
colei  water  upwellee!  mar  the  coast  act'  a'  dve.  tracing  features  in  the-  flow  field. 

I' he-  alia'  reinin'  a  much  different  impression  of  the-  upwelling  circulation  than 
that  obtained  from  the  mean  field  measurements  discussed  in  (  -hapt e-r  (>.  T  he  mean  fieTd 
showed  offshore  flow  ne  e  urring  in  a  surface-  Lkrnati  layer  at  velocities  of  order  10  cm  sec 
or  less.  I’lie  svnoplie  maps,  however,  reve’a)  greatly  enhanced  offshore-  flow  occurring  in 
narrow  band-  along  the-  eoa't.  Also  observed  in  drifter  records  from  the  same  period. 
Dav  i-  llflstal  Ini'  dubbed  I  lie-se  feature-s  'Lcpiirt'".  On  27  April  19hl  (fig.  7.1).  a  15  km 
wide  sipnri  i-  seen  centered  25  km  south  of  Pt .  Arena,  ('old  water  dyes  a  portion  of 
iliis  sepnrt .  but  both  the  l)Al.  thermistor  record  and  the  satellite  image  show  that  the 
sepiirt  i'  somewhat  w  ider  than  the  cold  tongue.  Compensating  onshore  flow  at  depth 
doe-'  not  occur  within  the-  DAI  range  of  15Um;  this  flow  is  truly  three-dimensional.  To 
the  north  arul  south  of  the  sepiirt.  warm  surface-  wate-r  accompanies  shoreward  flow. 
Another  sepiirt  carrying  upwelled  water  seawarei  crosses  the  Ross  line  in  the-  south  (cf. 
I  ig.  (>.2al  Despite  generally  strong  wind  forc  ing,  alongshore  currents  are  weak,  and  the 
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<  on-la'  jet.  winch  dominates  l  Ik  mean  flow,  i-  absent. 

A  week  later,  tin  region  wa-  resun  ey  ed  (fig.  7.2)  during  a  period  of  strong 
wind  ft  >rt  ing.  A  much  more  classical  flow  i'  observed  in  t his  survey.  with  currents  polar¬ 
ized  alongshore  over  most  of  t  Ik  shell.  However,  off  I’t.  Arena  tin  currents  turn  sea¬ 
ward.  again  earning  upwelled  water  away  from  the  coast.  Surface  drifters  deployed  in  a 
lint  off  I’t.  Arena  dramatically  echo  1 1 1  i  ^  observation  (set-  fig.  a  in  Davis  (19Ma)|. 
Nearly  Kin  km  from  shore,  st  rtnig  narrow  current'  can  be  seen  deforming  tin  tempera- 
tun  field  in  tin  southern  portion  of  the  survev. 

I  ig  makes  it  clear  that  energetic  eddies  centered  over  the  continental  slope 
tan  reach  onto  tin  shelf  and  deform  the  temperature  field,  sweeping  upwelled  water  from 
tin  shelf  and  even  recirculating  it  back  again.  This  survev  appears  to  show  a  counter¬ 
rotating  eddy  pair,  joined  along  I’t.  Arena.  The  strong  on  offshore  currents  from  these 
edd  ies  wen  still  apparent  when  the  area  wa'  resurveyed  (fig.  7.1 ) . 

fig'.  7.7  through  7.11  chronicle  tin  development  of  the  temperature  and  current 
fields  over  a  two  week  period  in  mid-July  19*2.  The  winds  varied  greatly  at 
MMU)  KiOl.'i  during  thi-  period  (fig.  (>.l)  and  shipboard  winds  (Buyer.  ft  al  19><1)  show 
strong  variability  during  the  course  of  each  survey,  fven  so.  strong  features  un¬ 
observed  in  the  currents,  measured  over  several  days,  which  correspond  well  to  features 
in  the  sea  surface  tempera!  urt.  measured  in  minutes  by  the  satellite.  Two  cold  tongues 
arc  observed  in  the  images  of  9-K>  July  IfWi.  The  northern  longue  is  swept  out  to  sea 
along  the  southern  edge  of  a  strongly  sheared  current  feature  which  hugs  the  coast  north 
of  I’t.  Arena.  The  southern  tongue  exhibits  weaker  flow  along  it  axis  and  weaker  shear 
across  it.  By  the  last  survey  (fig.  7.11).  the  southern  tongue,  and  the  currents  associ¬ 
ated  with  it.  have  disappeared.  The  development  of  a  cyclonic  eddy  between  the  two 
cold  tongues  can  be  seen  in  the  IR  and  also  in  the  current  measurements  of  f  igs.  7. Ill 
and  ill.  Ibis  eddy,  (entered  over  the  continental  slope,  dominates  the  dynamic  topog¬ 
raphy  during  the  lb-19. Inly  survey  (Buyer,  ft  al.  19M).  and  one  surface  drifter  was 
trapped  by  the  eddy  for  nearly  h  days,  completing  .'{  circuits  around  its  core. 

hi  summary,  the  maps  show  a  current  field  which  varies  strongly  in  both  x  and 
y.  and  which  only  rarely  resembles  the  smooth  average  field  described  in  Section  (>.2. 
Instead  ol  occurring  in  a  simple  surface  f.kman  layer,  significant  offshore  transport  of 
cold  coastal  water  i'  seen  iri  very  active  regions  of  short  horizontal  scale.  A  jet  in  the 
alongshore  current,  the  most  prominent  feature  seen  in  the  mean  field,  is  seen  only  occa¬ 
sionally  in  the  individual  maps.  Strong  current  fluctuations  from  the  mean  Held  are 
present  at  a  wide  range  of  scales,  and  coherent  eddy-like  structures  centered  beyond  the 
shelf  break  carry  water  upwelled  at  the  coast  far  out  to  sea. 


I  wo  striking  leal  lire.'-  of  I  lie  DAl.  map'-  (Icmtw-  eli-e Mission. 


7.4.  Wind  Relaxation  of  April  1982 

Hie  s  1 1  r  v « •  \  from  Leg  (>  of  ('Ol)l.  2  affords  a  well  defined  example  of  a  relaxation 
of  tin  current  field  on  cessation  of  strong  upwellmg  favorable  wind  forcing,  following 
sustained  ecpiatorward  winds  which  commenced  on  1-1  April  19*2  and  apparently  trig¬ 
gered  the  spring  transition  (Hover,  .''obey  and  Smith  1979)  to  the  upwelling  dynamical 
regime  in  ih<  (  '  ( >1)1.  region  (Lent/  and  Winant.  in  preparation),  tin  winds  over  the 
entire  (‘<>1)1  2  array  calmed  dramatically  <>n  19  April,  arid  remained  calm  until  27 

\|>rii.  I  In  coastal  current  field  was  mapped  by  shipboard  DAI.  during  29-2-1  April. 

I  In  satellile  image  from  this  period  (f  igs.  7.7  and  7.(>  in  the  atlas)  shows  a  wedge  of 
warm  surface  water  close  to  the  coast,  with  cooler  water  King  offshore.  Simultaneous 
Hydrographic  measurements  ( I  leise  hbein.  (Hilbert  and  Iluycr  19*2)  show  that  the  warm 
surface  water  L  accompanied  very  close  to  shore  bv  a  thin  layer  of  low  salinitv.  The 
Hiis'ian  Kiver.  flowing  at  an  average  rate-  of  107  m'se-c  1  during  this  period  (Markham 
>t  it!..  19*1 1  is  a  possible  source  for  the  fre-sh  water,  but  is  probably  too  small  to  explain 
i  in  extent  of  tfn  anomalously  warm  water. 

the-  observed  flow  field  L  remarkable  (figs.  7.7  and  7.(>).  Near  the  coast.,  strong 
poievvard  flow  is  associated  with  the  band  of  warm  water.  This  flow  occurs  in  the 
absence  of  any  significant  wind  forcing,  further  offshore,  a  strong  equatorward  jet  is 
present  f  nlike  the-  other  reali/.at  ions  of  the  flow  field  shown  in  the  atlas,  this  flow  is 
strongly  polarized  in  the  alongshore  direction.  The  depth  structure  of  the  alongshore 
flow  across  a  series  of  ('Tl)  transects  covering  a  2(>  hour  period  is  shown  in  Fig  7.12. 
displayed  from  north  to  south  Nearshore  poleward  flow  is  observed  in  each  transect, 
albeit  strongly  reduced  at  the-  Irish  (Hulch  Line  north  of  I’t.  Arena.  The  observed 
cyclonic  horizontal  shear  L  strongest  along  the-  Arena  Line-,  where  dv  dx  is  somewhat 
large  r  than  the  (  oriolis  parameter  f.  The  observed  vertical  shear  structure  is  interesting. 
I  hi  poleward  flow  L  relatively  barotropic.  while  the-  offshore  zone-  shows  strong  vertical 
shear  o!  19  ‘sec  1  and  more.  Ibis  suggests  that  the-  nearshore  zone  is  barot  Topically 
forced,  as  bv  an  alongshore  sea  level  gradient.  The  barotropic  component  of  the  flow 
can  not  be  determined  hv  drograplnc  ally  .  As  a  result,  the  surface  dynamic  height  maps 

I  roil  i  this  period  fail  to  delect  the  poleward  flow  near  the  coast  ( lleischbein.  (Filbert  and 

II  uyer  19*7). 

I’he  persistence  of  the  strong  structures  discussed  here  is  also  noteworthy.  Dur¬ 
ing  thr  cruise,  the  Central  line  was  surveyed  1  times.  Fig  7.12.  which  shows  offshore 
profiles  of  the  alongshore  current  at  2*m  depth  for  each  survey,  shows  that  the  cross- 
shelf  profile  of  v  changed  very  little  in  more  than  2  days  within  20  krn.  of  the  coast; 
hey  ond  27  krn  offshore,  the  ecpiatorward  jet  in  \  apparently  sharpened  somewhat 


betvvccii  liu  >('(<111(1  and  third  survevs.  ''(•(lions  from  cacti  survey .  1  lg  7.1-1.  an  a  No  reia- 
t  i \  * •  i >  unchanged  over  .'i  (lavs,  except  for  some  weakenin';  of  t  tie  field  >een  in  I  tie  last  sur- 
'  n  Any  theory  which  seeks  in  account  for  these  relaxation  event*,  must  explain  thi- 
long  decay  I  inn  . 

All  interest  hi"  footnote  to  thi>  discussion  coin  erii-  the  near-surface  shear  field 
during  relaxation  event-  (I  ig  7 . 1  r»  1 .  Prior  to  tin  spring  transition,  the  near  surface  shear 
field  show  *  strong  amplitude  fluctuation'  at  diurnal  and  higher  frequencies.  for  most  of 
tin  post -t  ran>it  ion  data,  these  fluctuation'  an  si  rough  damped.  During  periods  of  wind 
relaxation,  however,  tin  u-ualh  well  mixed  surfaci  laver  of  water  re-t  rat  ifie-.  allowing 
near  -urfaii  'hear-  to  he  supported,  and  tin  ampin  ud<  of  tin  shears  increases  drama!  i- 
call>.  for  ('.‘I.  tin  rest  ratification  occurs  not  only  in  temperature  (fig  7.1’>|  but  also  in 
-ability.  a-  can  be  seen  from  tin  h>  <1  rographii  data  ( f  leisi  hbein.  (iilbert  and 
Miner  lfx.'.i,  So.  (or  thi-  site,  advection  possibly  supplemented  by  local  surface  heating 
lorin  a  strong  surface  signature  iti  tin  density  field  during  relaxation  events.  When  this 
(tensity  structure  is  present,  strong  vertical  shears  can  be  supported,  and  are  observed, 
near  tin  surface.  Strong  near->u  rfa<  t  vertical  shear  is  usually  associated  with  strong 
wind  forcing,  not  with  its  ab-cnce.  Tin  energy  is  concentrated  in  motions  with  diurnal 
am:  semi-diurnal  period',  suggesting  tidal  forcing. 

f.xaminatioti  of  feature*  -eeti  during  relaxatioti  events  has  provided  insight  into 
•  i'  pert  s  of  t  in  momentum  baiarnc.  It.  tin  next  section,  we  shall  examine  features  which 
may  have  importance  to  tin  wav  s  in  which  nias-.  heat,  salt  and  nutrients  are  exchanged 
between  coastal  and  offshore  water' 

7.5.  Jet  Surveys,  July  1981  and  July  1982 

Tin-  features  in  tin  (‘01)1.  measurements  which  have  drawn  the  most  interest 
art  tin  prominent  tongues  of  coid  water,  a  part icularlv  active  example  of  which  is  shown 
iti  1  ig.  7.1b.  first  observed  in  satellite  imagery  (Bernstein,  ct  al.  1972.  Breaker  and  Cil- 
li. uni  19x1.  I  ragan/a  ft  al.  I9M.  hcllv  19x.‘i].  these  features  ran  extend  from  the  coastal 
iipwrlling  /one  to  several  hundred  kilometers  offshore.  In  the  example  shown,  at  least  t 
-trong  (old-water  baud'  art  *ccii  in  t  h<  *00  km  stretch  of  coastline  from  Cape  Merido- 
(  mo  to  ft  Conception.  The  potential  importance  of  these  structures  to  the  cross-shore 
mass  and  heat  haluine-  ha'  been  noted  by  Dav  is  (198-lb).  If  these  modes  of  proper!  v 
transport  an  significant,  their  Iran-lent  nature  and  strong  spatial  variability  imply  that 
an  iinder-t and ing  of  sin  h  balances  based  on  moored  data  alone  will  be  difficult,  at  best. 

I  In  cariie't  finding'  in  the  fir-t  week-  of  ('OI)K  with  both  DAI.  (figs.  7.I-7.-1) 
<ind  surface  drifter-  (Davis.  I9x.‘i)  were  that  unexpect  edlv  large  cross-shelf  currents  on 
narrow  scale-  were  associated  with  tin  cold  water  /ones  observed  in  the  satellite  images 
of  Kelly  (  19x21.  \*  a  re-nlt.  w  <  undertook  DAL  and  hydrographic  surveys  of  the  cold 


witter  le-atures  in  .1  ■  1 1  v  o(  19s)  and  again  in  ,lnl\  I9SJ.  (  <  >  r  1 1  r  1 1  >  1 1 1  inn-  in  tin  I9S1  surre> 
were  made  hr  .lam  llurer  ami  Martin  Oiirera  o’’  (>-l  (hydrography;.  ami  supplemcn- 
iar>  information  useful  in  cruise  planning  .iini  laier  anah-m  weri  provider!  hr  hathrrn 
kcllr  <>!  S|()  tsaicliite  !  rt  iii”<  ami  analrsi-  and  Larry  Breaker  of  \L>>  (satellite 
analr  -i-  ; .  ’!  in  19'''.’  Miner  drew  on  tm  re-soure  e-  of  ,1am  timer  I  nr urography  i.  am: 

I’ierre  |- lament  .  I.ihe  Washburn  ami  l.arrr  Arini  of  SIC)  lunelcrway  t  he  rmosalinograph . 
tow-ro  ('ll)  measurement-  ami  -hon  to  ship  t  ransmtssioi,  of  satellite  image- i. 
Alt  hough  ttu  latter  group  -  primurr  interest  was  Hi  small  scale  mixing  proe  esse-  at  tin 
tioimdarie-  of  t  he-i  (old  water  structure—..  their  satellite  image-  and  underwar  dal;-, 
allowed  much  more  ifleetiri  planning  of  survey  strategies  for.  and  scientific  understand¬ 
ing  of.  tin  mcso-rale  -t  rue  tun—  w  im  li  were  tin  primary  interest  of  tin  DAI.  and  hreiro- 
gr.ipli ic  group- . 

\1  ap-  of  tin  currents  ai  '."'in  depth  from  tin  two  surveys,  tin  first  ((inducted  ot; 
1- 19  .1 11 1 \  !9sl.  ihi  second  on  22-2F  ,liil\  I9s2.  arc  shown  in  Figs.  7.17  and  7  lv.  Both 
-urrer-  show  a  tongm  of  cold  wider  riding  <>ti  \  err  swift  (up  to  7a  cm  sec  i  offshore 
(iirrent-  from  their  apnaren!  point-  of  origin  in  tin  etiastal  upweliing  /one.  exte-neiing 
orer  'Jan  Km  out  to  s(  a.  \\  itliin  tin  itroad  i  (H7~>  km  wide)  tongues  of  cold  water,  a  nar¬ 
row  i  <  >  ; ’>  km  wide  hand  of  crei  i  oider  water  lies  toward  the  northern  edge  Some  er  i- 
denee  of  recirculation  hail,  toward  tin  (  oa-l  in  the  southern  diffuse  portion  of  t  he  longue 
e  seen  ii,  i  he  1 9'- 1  surrer  1  iamen:  (personal  communication)  has  noted  thi-  recircula¬ 
tion  in  seepiene  c-  of  satellite  image-  for  the  19VJ  surrer;  unfurl  u  nat  elr  in  1 9  s  J  our  in 
sih;  measurements  did  not  i  ricnil  liir  enough  south  to  sample  this  region. 

The  surrer  of  JJ-Jh  .luir  I9s2  was  i  h<  more  successful  at  obtaining  sections  of 
oceanic  variability  across  the'  cold  tongue,  primarily  because  of  tlie  shipboard  availability 
ol  sal  ell  it  e  image-  from  a-  late  as  JJ  July  J.'iOO  l  T.  By  the  time  the'  surrer  was  eon- 
dinted.  cloud-  had  obscured  tin  area,  and  the  detailed  structure  of  the  field  had  altered, 
but  the  data  Irorn  tin  surface  t  hermosalinograph  allowed  the  actual  locations  of  featurc- 
il  t  hi  sorr  cr  to  be  blent  ified. 

flu  ship  track  lor  t  lie  I9V2  surver.  annotated  with  t  lie  dar  of  the  rear  and  time 
I  I  ).  t-  show  i:  it  I  tg  7  19  (a  gap  in  the  DAI  data  occ  urred  on  21  July  (day  205)  from 
I  I A  0  -  191  I  l  I  i  some  ele  eon  crossings  were-  made  of  t  lie  thin  band  ol  coldest  water  at 
tm  northern  edge  ol  the  cold  longue.  When  the  underwar  t  hermosalinograph  data  from 
tins  period  i-  plotted  as  a  T-s  diagram  (I  ig.  7.20).  a  strong  separation  of  1 -S  behavior 
a  i  ri  is-  the  told  tongue  is  found  The-  left  branch  of  the  curr  e.  made  up  of  samples  within 
and  north  of  the  coldest  ribbon  of  water  represents  relative!)  fresh  water  whose  saltnitr 
raru-s  lit  tic  orer  a  temperature  range  of  nearly  •’>  <  .  I  he  right  hand  branch  has  higher 
salinity  and  lie-  south  of  the-  coldest  ribbon.  From  ('ll)  sampl  on  stations  in  the 

I9s|  surrer.  Olirera  (I9s2)  a|sn  postulate-  a  separation  of  surl.u  e  water  types,  and 
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attribute-  tin  northern.  Ire-her  water  ma~- 
plunn 


dilut  cd  remnant  of  (  oiumlua  K ;  \  < 


It  order  in  rr-nlu  I  lit  I )  \  I  ciirrenl-  min  ><■(  lime  -hovving  their  -triieiurc  win. 
depth  and  |  >e  *—  1 1  i 1 1 n .  a  (  oortli nat e  -v-icm  must  be  chosen.  I  Ik  i  limn  i-  no:  obviou-. 
However,  -mce  i  Ik  ie:  t  hangc-  it-  orient  aliori  al  eat  i.  (  ro--ing  !  (win  t  adopt  ed 
m  rie'iei  aiongje:  coordinate-  l\.\,  with  current'  lu.vl  defined  to  in  across  alone  tin 
direction  o:  ma--  t  ran-port  me-a-uri-d  w  it  n  t  tic  I )  \  I  aero--  each  transect.  Section-  were 
:  lien  pint  ted  a-  tunct  ton-  o!  t  lit  ae  ros-ie;  eoordinati  \  Because  t  (it  jet  alvvav  -  point  - 
ofl-iiort  .  the  northern,  low  -alinitx  water  i-alwa\-  at  lush  \:  tin  origin  of  \  i-  arhit  rarv . 

1  or  eac  h  set  t  ion.  current-  were  averaged  over  Inn-  of  a  km  in  x.  I  he  aiongjet  current' 
nave  heel  contoured  at  Id  cm  set  intervals:  the  acrossiet  current-  art  show  i.  a-  arrows 
u-mg  a  -tali  of  >  kin  in  \  for  .'idem  sec  in  u.  Shown  above  each  section  are  tin-  con¬ 
current  I.  >.  and  den-itv  (a-  a  tit  a  m  depth  from  the  underwav  t  hermosaiinograph . 

•\t  each  e  ro--ittg.  a  surface  intensified  jet  in  v  is  observed,  with  peak  current- 
ranging  iron,  over  TO  cm  set  to  less  than  -Id  cm  sec.  Within  tin  jet.  tin  aiongjet  flow  is 
20  (  in  -el  or  mon  to  at  leas;  Idd  in  depth.  I  In  iet  appear-  to  acliievt  its  maximum 

veiocitv  not  a!  I  In  -harp  t  emperat  u  rt  front  when  tin  (oldest  water  i-  found,  but  onlv 

alter  tin  temperature  bee  in-  t , ,  rise  again  at  larger  \.  i.  e..  within  the  northern  water 
ma'.-.  Tin  data  show  that  tin  dramatic  temperature  front  seen  in  tin  satellite  images  is 
-tronglv  -alinitv  compensated,  so  that  surface  den  sin  change-  lilt  it  across  it.  while  tin 
slow  rist  it,  surface  temperature  north  of  the  front  t-  ac<  ornpaniee!  I  •  v  fair!)  constant 
salmitv.  rt-ulting  in  a  gentle-  densiiv  gradient  This  den-itv  gradient  will  in  turn  he 

•it  e omoanied  In  a  rise  in  geosirophie  alongjet  current,  Thu-,  at  ica-t  near  the  surface, 

tin  saliniiv  structure  contributes  independent!)  to  the  densiiv  Held,  and  estimates  of 
geostrophn  flow  based  on  temperature  measurements  alone  will  hi  in  i  rror. 

Within  tin  l)  \l.  field  of  view,  tin  measured  aiongjet  transport  varies  from  tl.i 
to  2. .‘I  vv  (  I  sv  Id’  Hi'  see  ),  with  an  average  of  l.a  Sv  over  eompiete  transects.  I  hi- 
wih  he  an  unde-rest  imatc  of  tin  ma--  transport  stnec.  in  man)  ease-,  tlu-  measurement - 
did  not  reach  deep  enough  nor  far  enough  north  to  define  tin-  boundaries  of  tin  flow. 
Nonet  hi  ie--.  thi-  represents  a  stunning  amount  of  water  flowing  offshore,  for  com- 
parison.  tin  offshore  l.kman  transport  resulting  from  a  surface  windstri-s-  r  of 
1  i  i  v  im  <  iii'  i-  r  /if  I  m“  see.  |f  j  B  i  ~  transport  were  to  he  provided  bv  a  cum  ergenc ( 
of  l.kttirit.  transport  Irotti  the  coa-tal  zone,  it  would  require  gathering  all  the  l.kman 
transport  generated  bv  a  la  elv  tie  cm*  wind  over  1(100  km  of  coastline.  Clear!)  the  hulk 
of  the  water  travelling  offshore  i-  not  being  provided  bv  sirupit  convergence  of  water 
u  pw  el  led  at  the  coast.  Moreover,  anv  flow  mitiallv  restricted  to  the  shelf  would  have  to 
be  prodigious  indeed  to  suppiv  this  much  mass;  the  total  diversion  of  a  current  averaging 
7.')  cm  see  over  the  entire  breadth  and  depth  of  a  20  km  wide  shelf  1  OOrti  deep  would  be 
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Fig  7.-1  :  Ship  surve\  from  23  May  0100  to  25  May  0600.  1981.  NO  A  A  6  image*  iron 
20  Ma>  1550  (same  as  Fig.  7.3). 


Fig.  7.5  :  >hip  survey  from  20  April  1600  to  22  April  1200.  1982.  NOAA7  image 
from  24  April  2210  (note  delay  from  survey  time).  This  flow  occurred  during  a  cessation 
of  winds  following  a  period  of  strong  southward  wind  forcing,  and  is  discussed  in 
Section  7.3. 


i 


Fig.  7.C  :  Ship  survev  from  22  April  1401)  lo  24  April  1130.  1982.  NOAA7  image 
from  24  April  2210  (same  as  Fig  7.5). 


Kig.  7. 8  •  Ship  survey  from  11  July  1100  to  Ki  July  1400.  1982.  NO  A  A"  image  from 
9  July  2211  (same  as  Fig.  7.7) 


Fig.  7.9  :  Ship  survey  from  14  July  1400  to  16  July  0920.  1982.  NOAA7  image  from 
14  July  2254.  Wind  forcing  was  very  strong  and  increasing  during  the  survey  period. 
Note  the  strong  current  magnitudes  and  shears  across  the  northern  transects. 


Fig.  7. JO  :  Ship  survey  from  16  July  1730  to  18  July  0845.  1982.  NOAA7  image  from 
15  July  2241.  Survey  was  performed  from  south  to  north:  winds  were  strong  to 
moderate  along  the  Central  and  North  lines,  but  had  calmed  to  airs  during  the  surveys 
of  the  Arena  and  Irish  Gulch  lines.  Northward  flow  at  the  inshore  ends  of  the  CTD 
lines  is  seen,  although  the  flow  at  the  offshore  end  of  the  two  northern  lines  remains 
strong  and  strongly  sheared.  Note  the  developing  eddy  along  the  North  line. 
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Fig.  7.11  :  Ship  survey  from  19  July  1900  to  22  July  0400.  1982.  NOAA7  image  from 
22  July  2259.  Winds  were  calm  during  the  survey  of  southern  3  lines,  but  were 
moderate  to  strong  along  the  northern  3  lines.  The  eddy  along  the  North  line  has 
developed  strongly.  Davis  (1983)  shows  tracks  of  surface  drifters  w'hich  circulated  in  this 
eddy  for  a  week  and  more,  making  several  circuits. 


Appendix 


DAL  SAMPLING  AT  JIYDKOORAFMIIO  STATION'S 

This  appendix  shows  tin  times  (day  of  tin  year)  at  which  DAL  data  were  col 
letted  at  each  hydrographic  station  shown  in  Fig.  (>.l.  Occupations  separated  by  les 
t  han  I  day  are  not  shown. 


T-S  diagram  for  jet  survey,  204/08:15  to  206/19:30 


Fig.  7.20  :  T-S  diagram  from  underway  thermosalinograph  during  jet  survey  (data 
kindly  made  available  by  Pierre  Flament.  Libe  W  ashburn  and  Larry  Armi).  Sampling 
depth  is  approximately  5m. 


Fig.  7.19  :  Ship  track  during  1982  offshore  surve\.  indicating  locations  of  sections 
through  the  cold  water  tongue.  Clouds  obscured  the  area  during  the  actual  survey  time, 
and  the  location  of  the  front  and  detailed  structure  of  the  temperature  field  were  not 
available  from  satellite  III. 
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Fig.  7. IS  :  Time  series  of  moored  measurements  showing  alongshore  component  of 
wind  stress  r  at  NDHO  "10(113,  ami  temperature  difference  and  Ov/ftz  between 
instruments  at  .r>  and  10m  depth  at  ('3  vs  day  of  1982. 
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occupations  vs.  day  of  year. 
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